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THE RECEPTIVITY OF BOUNDARY LAYERS ON BLUNT BODIES 
TO OSCILLATIONS IN THE FREE STREAM 


By 

Chester E. Grosch* 
INTRODUCTION 


The prediction of the position of transition on a body is one of the 
most. difficult problems in fluid dynamics. Linear stability theory with the 
"e n " criterion (see Mack, 1977, for an extended discussion) is probably 
the most widely used method of predicting transition. The major weakness of 
this method is that it employes a relative measure of flow instability; the 
growth of some measure of the disturbance size, whether amplitude, energy, 
or whatever, relative to an unknown initial size of the disturbance. 

It is clear, as pointed out by Mack (1977), that the position of tran- 
sition must depend on the absolute size of the disturbance, and thus on the 
initial size of the disturbance as well as the growth rate of the unstable 
disturbance. If there were no disturbance at all, that is the initial size 
were zero, the stability or instability of the flow would be irrelevant. 

The flow would remain laminar, and transition would not occur anywhere on 
the body unless, of course, the flow separated. On the other hand, if the 
initial disturbance were very large, say of the order of the mean flow 
speed, then transition would, presumably, occur very near the front of the 
body, perhaps in the immediate vicinity of the stagnation point or line. 

A rough estimate of the magnitude of the initial amplitude of the 
Tollmien-Schlichting waves at the beginning of the region of instability can 
be made if it is assumed that the "e 11 " criterion is an approximately valid 
transition method and if one accepts the estimate of Klebanoff, Tidstron, 
and Sargent (1962) that, at the beginning of transition, the rms value of 
the perturbation velocity is about 1.5% of the free stream speed. Since the 

*Professor, Department of Oceanography, Old Dominion University, Norfolk, 
Virginia 23508. 



initial proposal of the "e n " criterion by Smith (1956) and Van Ingen 
(1956), extensive studies of transition data have shown that n varies from 
about 2.5 to 10.0, depending on the test facility. Taking the disturbance 
amplitude at transition to be 0.015 of the free stream speed, and n to 
vary from 2.5 to 10.0, it follows that the amplitude of the disturbance at 
the beginning of the region of flow instability is in the range 10“^ to 10“^ 
of the free stream speed. This suggests that a very weak coupling of free 
stream disturbances to To 1 lmien-Schl icht ing waves in the boundary layer may 
be extremely important in the initiation of the transition process. If the 
free stream disturbances have amplitudes of 0(e), even disturbances in the 
boundary layer of 0(e 2 ) may play a role in initiating transition. 

In order to incorporate the influence of the flow environment in sta- 
bility prediction methods, and thereby provide a rational basis for these 
methods, it is necessary to be able to calculate the initial size of the 
disburbance caused by external, i.e. free stream, disturbances. There are 
four classes of free stream disturbances which may be responsible for the 
generation of the Tollmien-Schlichting waves which lead to transition 
(Obremski, Morkovin, and Landhal, 1969). They are: mean flow unsteadiness, 

free stream vorticity, sound, and entropy fluctuations. If the fluid were 
truly incompressible, only the unsteadiness of the mean flow and free stream 
vorticity could exist and initiate the transition process. In reality, 
sound waves are always present. Entropy fluctuations appear to occur only 
at high supersonic speeds. In a recent study, Harvey and Bobbitt (1981) 
have examined the anomolies between wind tunnel and flight transition mea- 
surements. They state that the experimental results show"... the dominance 
of free-stream disturbance level on the transition process from beginning to 
end." Further, they conclude that there is a question of"... whether sound, 
unsteadiness, or spectral peaks is the most influential on boundary layer 
receptivity,..." Morkovin (1978) has argued that, in medium to low subsonic 
flows, sound waves can be modeled quite accurately by a time varying 
oscillation superimposed upon a steady mean flow, provided that the 
wavelength of the sound wave is much larger than the characteristic length 
scale of the body. For a sound wave impinging on the nose of a body this 
characteristic scale length would appear to be the radius of curvature of 
the body at the forward stagnation point. 
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Recent studies, (Salwen, Kelly, and Grosch, 1980; Grosch and Salwen, 
1980; Salwen and Grosch, 1981) have been addressed to the problem of calcu- 
lating the amplitudes of the discrete and continuum eigenfunctions of the 
linearized, parallel flow stability problem. The results of these studies 
can be summarized as follows: first, the proper inner products for ortho- 

gonality of these eigenfunctions were found; second, it was shown that the 
set of discrete and continuum eigenfunctions for both the temporal and spa- 
tial problems is complete; and third, a procedure for calculating the ampli- 
tudes of the discrete and continuum eigenfunctions was found. In order to 
calculate these amplitudes, the disturbance must be known, as a function of 
time, on a plane perpendicular to the boundary. 

In order for this procedure to be applied it is necessary for the dis- 
turbance to be small so that linearization is valid; that the parallel flow 
approximation be valid, at least approximately; and that the form of the 
disturbance be known throughout the boundary layer and free stream at one 
location on a two-dimensional body. Even if the disturbances are small, 
this theory cannot be applied in the immediate vicinity of the leading edge 
of a body because the flow in that region is not even approximately 
parallel. A basic problem is to determine the flow within the boundary 
layer caused by disturbances propagating and being advected by the mean flow 
toward the front of a blunt body. 

The response of the boundary layer on a body to oscillations in the 
free stream has been the subject of a number of studies beginning in the 
early 1950's. Some representative examples are the work of Moore (1951), 
Lighthill (1954), Stuart (1955), Rott (1956), Lin (1956), Glauert (1956), 
Carrier and Di Prima (1957), Gibson (1957), Watson (1959), Lam and Rott 
(1960), Rott and Rosenzweig. (1960), Sarraa (1964), Stuart (1966), and 
Ackerberg and Phillips (1972). Riley (1975) and Schlichting (1979, Chapter 
15) give comprehensive reviews of recent work in unsteady boundary layer 
theory. None of these authors have considered the application of their 
results to the boundary layer receptivity problem. 

These studies have considered either the unsteady flow in the vicinity 
of a stagnation point or on a flat plate. The flat plate problem presents 
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great difficulties because of the singularity of the boundary layer equa- 
tions at the leading edge of the plate. Carrier and Di Prima (1957) studied 
this problem within the framework of a linearized, Oseen approximation, to 
first order in the amplitude of the oscillation. They found that, far from 
the leading edge of the plate, the time dependent part of the flow is a 
Stokes shear wave. Ackerberg and Phillips (1972) studied the same problem 
using asymptotic, again to first order in the amplitude of the oscillation, 
and numerical techniques. They also found that, far downstream of the lead- 
ing edge, the flow develops into a nearly Blasius steady flow and a Stokes 
shear layer imbedded within the steady boundary layer. The most surprising, 
and unexpected, result of Ackerberg and Phillips is that "... most flow 
quantities approach their asymptotic values far downstream through damped 
oscillations. " 

Grosch and Salwen have found a solution of the Navier- Stokes 
equations for an incompressible stagnation point flow whose magnitude 
oscillates in time about a constant, nonzero, value. (See Appen- 
dix B.) The analytic solutions, to third order in the amplitude, obtained 
using inner and outer expansions, are in complete agreement with the results 
o.f numerical integrations. The mean flow is the steady stagnation point 
flow plus second, and higher, order flows driven by the Reynolds stress of 
the oscillatory components. The oscillatory flow is that of the fundamental 
and all of the higher harmonics of the fundamental. 

The mean flow is, to lowest order, independent of the oscillating flow. 

The structure of the oscillatory components of the flow depends on the fre- 
quency of oscillation, not on the distance from the stagnation point. As 
was expected, the oscillatory flow develops into a Stokes shear layer at 
high frequencies. The results of Carrier and Di Prima and Ackerberg and 
Phillips that the Stokes shear layer develops far downstream may be due to 
the use of the boundary layer equations in a region in which they are 
invalid. 

In any case, the flat plate with zero thickness and a "sharp" leading 
edge is a nonphysical model; all physical bodies have a blunt leading edge, 
at least if a continuum approximation can be applied to the flow past the 
body . 
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In this paper we examine one aspect of the boundary layer receptivity 
problem; that of the flow induced in the boundary layer on the forward por- 
tion of a blunt body by unsteadiness in the mean flow. As was mentioned 
above, this is also an approximate solution for the disturbances generated 
in the boundary layer by a long wavelength sound wave impinging normally on 
the front of a body in steady motion. This work is an extension of our 
previous study of the flow in the neighborhood of a stagnation point (Grosch 
and Salwen, Appendix B) . 


FORMULATION 


Consider the flow past a blunt, two-dimensional body. The usual bound- 

r\j 

ary layer coordinates are used: x is the distance along the surface of the 
body from the forward stagnation point and y is the distance from the body 

/v 

along the normal to the body. A solution is to be found for x 2: 0 of the 
time-dependent, two-dimensional boundary layer equations 


3 u 


3v _ 





+ 


o. 



3 x 


3y 




3 u 


3 u 

3 u 

3U 

, „ 3U 


+ 

u 

+ v — = 


+ U — 

3 t 


3x 

3y 

3 t 

3 x 


+ v 


3 2 u 

3 

3y 2 


(1) 

(2) 


with the boundary conditions 


u(x,o , t) = 

v(x,o,t) = 0, 

(3) 

A/ 

u(x,y,t) -> 

U(x,t) as y + 00 . 

(4) 


Here, (u,v) are the components of the velocity in the x and y 
directions, U(x,t) is the component of the potential flow which is tan- 
gential to the body surface, and v is the kinematic viscosity. It will be 
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assumed that the potential flow is a temporally oscillating flow, super- 
imposed on a steady mean flow; specifically it is assumed that, 

U(x,t) = U Q H(x/^)[l + e cos tot], (5) 

where U 0 is the magnitude of the mean potential flow, 
l is an appropriate length scale of the body, 

£ is the dimensionless amplitude of the oscillating component, 
and, 

H is a dimensionless function, giving the variation along the body 
of the tangential component of the potential flow. 

It will be assumed that the body is symmetric and is at a zero angle of 
attack, so that H is of the form 

H(x/£) = I a 2n+l ( x /^ 2n+1 > (6) 

n=0 

with the {aj} dimensionless constants. Specific examples will be 
considered below. 


METHOD OF SOLUTION 
Define dimensionless variables 


ll 

Xj 

(7) 

n = (a^Re) 1/2 (y/£) 

(8) 

t = ai U t/£, 
1 o 

(9) 

parameters, a Reynolds number 


Re = U £/v 
o 

(10) 
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and a Strouhal number, a dimensionless frequency, 


a = loA/CaiU ) . 

o 

Defining a stream function such that 

u = 3iJ»/3y, v = -3 tJ)/3x, 
by 

= [ (U o *)/( ai Re) 1/2 ] <K5,n,T), 
equation (1) is satisfied identically and (2) is 


(ID 


(12a, b) 
(13) 


3 3 

— — + 

3n 3 


-i 

a i 


3<(> 3 2 <t> 
3? 3h 2 


3<j) 3 2 <j> 
3q 353 ti 


3 2 <J> 

3t3ti 


= eaH(5) 


sin ax 


-1 


- a x H(5)H'(C)[l + e cos ax]- 


(14) 


The boundary conditions, (3) and (4) are 


3d> 3d> 

_L = _L=0 at n = 0, 

3 5 


(15) 


and 


— ■+■ H(5)[ 1 + e cos ax] as ri + ». 
3n 


(16) 


The tangential component of potential velocity, H(5), is a power series in 
5. It is clear, from (14), that <j> must also be a power series in 5, 
with each coefficient a function of h , X, and the parameter E. Each of 
these coefficients will be expanded in a Fourier series in x, with each 
Fourier amplitude a function of n and the parameter e. Finally, each 
Fourier amplitude will be expanded in a power series in e, with each 
coefficient of the power series a function of q . These functions of h 
are then found to be solutions of ordinary differential equations. 
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First, assume that 


t«,n,T) - [[«■ * l) - « 0 J * 2n+1 S 2n+1 F 2 „ + i (n " T) - C17) 

n=0 

Substituting (6) and (17) into (14) and equating the coefficients of 
successive powers of £ on the left and right sides of equation (14), the 
partial differential equations for the {F.(h,x)} are found. F^ satisfies 


F^ilL-fSFiV 8ilL. P( T ;e , 0) - q(t ;£ ,o) , (18) 

3q 3 1 3n 2 \3n / 3x3n 


where 


P(t ;e ,o ) = ecr sin ax , 


(19a) 


and 


Q(x ;e ,o ) 


1 + -L + 2e cos ox + — cos 2ox 
2 2 


Next, define the operators 


(19b) 


3 3 F . 3 2 F . 3 F 3 F . 3 2 Fj 3 2 F . 

D(F . ,Fj ; j) = i + Fi 1 - (j + 1) + j F. - i 

3 3n 3 3n 2 3 ti 3n 3ri 2 3 3x3q 


( 20 ) 


and 


3 F . 3 F . 3 2 F . 

G(F. ,F. ) = — - — i - F. . 

13 3n 3n 3ri 2 


( 21 ) 
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Then, 


D(F 3 ,Fi ;3) = I P - Q, (22) 

4 

D(F 5 ,Fi ;5) = I P - (1 + - C 5 ) Q + 8 C 5 G(F 3 ,F 3 ), (23) 

6 2 

D(F 7 ,Fi ;7) “ - P - (1 + C 7 ) Q + 9 C 7 G(F 5 ,F 3 ) + 15 C 7 G(F 3 ,Fs), (24) 

8 

D(F g ,F, ;9) = — P - (1 + C + - C ) Q + — C [48 G(F 7 ,F 3 ) 

10 9,1 2 9»2 5 9,1 

+ 1 12G(F 3 ,F 7 ) ] + 18 C G(F,- , Fj. ) , (25) 

9,2 D 3 

D(F n ,Fi ; 11) - 2_ P - (1 + C + C ) Q + C [30G(F 3 ,Fg) 

12 11,1 11,2 H,1 

+ 10G(F ,F )] + C [ 28G(F,- ,F 7 ) + 20G(F 7 , F,- ) ] , (26) 

9 3 11,2 b S 

with 

C 5 = a£/(a ia5 ), (27 a) 

Cy = 3385 /( 3 ! ay), (27b) 

C 9 , 1 = a 3 a? /(aiag) , Cg^ = a §/( a l^), (27c, d) 

c ll,l = a 3 ag /( a l a ll > » c ll,2 = a 5 a 7 /(a 1 ai L ) . (27e,f) 
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The boundary conditions for the ^2n+l^ are 


W 0 - T) 


= (3 f /an ) 
v 2n+r n=o 


= o 


( 28 ) 


and 


[ 2(n + 1) - 6 ] (3 F an) + 1 + e cos ot, 

L on JV 2n+l 


(29) 


as . n °° . 

We have previously found the solution for Fi up to, and including, 

terms in e 3 (Grosch and Salwen, Appendix B) . Although the 

boundary layer equations are being used in this study, it should be noted 

that 5 Fj (p ,t ) is a solution of the Navier- Stokes equations for oscillating 

stagnation point flow. It is easily seen that the form of the solution for 

each of the F„ , is the same as that for F, . Therefore, we have for 
2n+l 1 

j = 1,3,. ...11 


F . (n ,t ) = l 


2m 


„-0(n> * I £ 2nrt-k[ f 


n=0 


j ,2m. 


k=l 


j ,2m,k 


(n) cos (kcrr) 


+ g. „ , (n) sin (kOT )] 

J, 2m, k 


(30) 


It is quite straightforward, but somewhat tedious, to obtain the ordi- 
nary differential equations and boundary conditions for the f*s and g' s. 
These differential equations are given in Appendix A. 

One of the major problems associated with series solutions of this type 
is that of determining the radius of convergence. On the basis of the 
results from the study of oscillating stagnation point flow 
(Appendix B) , it can be concluded that the power series in e is rapidly 
convergent for e < 1, except perhaps at very low frequencies. If the fre- 
quency of oscillation is low, then a quasisteady approximation is quite 
accurate. Because of the rapid convergence of the series in e, for e 
small, only the terms through e 2 have been calculated. 
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A much more difficult problem is estimating the radius of convergence 
of the power series in £. It appears that there are no general results for 
determining the radius of convergence of the steady Blasius series. Van 
Dyke (1957) has shown that the steady Blasius series for a parabolic 
cylinder diverges for x/r Q > 0.62, where r Q is the radius of curvature of 
the parabola at the nose. The results given by Schlichting (1979) for the 
steady Blasius series for a sphere of radius r^ suggest that this series 

is convergent up to at least x / r Q ** tt/ 2. It seems reasonable to assume that 
the time dependent series solution given here is convergent for £ = x/l 1. 

RESULTS 


Numerical integration of the equations given in Appendix A yield the 

{ f. _ . ig. n Substituting (30) and (17) into (13) and using (12) 

j,2n,k j , 2nk 

gives, for the velocity components (u,v) , 


u/U = 


1 ' 2( " 

n=0 v 


2n+l 00 


+ 1) - « ] a. .£ 

on 2n+l 


X [e 2n £' 


m=0 


2n+l ,2m,0 


+ l £ 2m+k {q n+lj2mjk cos(kax) + g ; Ti ^ 1 7m v s in(kaT)}][, (31) 


k=l 


v/U - - (a- 
o 


- 1/2 “ ( 

1 Re) l <[2(n 
n=0 { 


2n+l ,2m, k 


+ 1) - 6 ](2n + 1) a ,£ 

on 2n+l 


2n 


2ntf-k 


l [ £ , o « + l {f 0 , „ , cos (kOT) + 

n L 2n+l,2m,0 . . 2n+l,2m,k 

n=0 k=l 


+ g 2n + l,2m,k sin (kIT)}] ‘ 


(32) 


The velocity field depends on the body shape, the {aj} , the 
Reynolds number. Re, and the dimensionless frequency, o. It should be 
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noted that the effective Reynolds number and frequency depend, to some 
extent, on the body shape because the velocity scale is a l^ Q > see 

equations (8), (10), (11), and (13). It is clear from (31) that, as usual, 
the tangential component of the velocity, u, has no explicit dependence on 
the Reynolds number which only affects the scaling in the normal (h) 
direction. Again, as usual, the normal component, v, is proportional to 
Re - !/^ addition to the scaled dependence of the normal coordinate. 

. o 

Through the terms in e , 


u/U = u (5 ,ti ) 
o o.o 


2 

+ e u 


2,o 


(5 ,n) + e[u (1 1 ) (5,n) 

o,l 


cos OT + 


(2),_ x 
u (§ ,n) 
0,1 


• 1 2 r (1) 0 (2) . i 

x sm ax + e I u _ cos 2 ot + u „ sm 2arl , 
o,z 0,2 


(3 3) 


v/U = - (a Re)~ 1/2 { v (5 ,h ) + e v (e ,r \ ) + e[ v^ ^ (£ ) cos ot 

o 1 o,o 2,o L o , 1 


+ v^ 2 j sin or] +e [v^(?,ri) cos 2ot + v^ 2 ^(5,n) 
o,l J L o,2 o,2 


(34) 


sin 2 jt]} 


with , 


u. « 1 [2(n + 1) - 6 ] a. ? 2n+1 f; . . (h), 

i,o „ on J 2n+l 2n+l,i,o 

j n =0 J 


(35 a) 


u (U - I [2<n * 1) - « J C 2 “ +1 C 




n=0 


on 2n+l 2n+l,o,j 


.(n). 


(35b) 


u (2 ? “ I [ 2(n + 1) - 6 ] a ? 2n+1 g' . (n), 

o,i _ L 0 n J 2n+l 2n+l,o,i 

J n=0 


(35c) 
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» 


- - I [ 2(n + 1) - Sj (2n ♦ 1) a S 2 "^, j ,„(» > . 

J n=0 J 


(35 d) 


v Cl ? - X [ 2(n + 1) - 6 1 ( 2n + 1) a„ . , £ 2n f _ , (n ) . 




n=0 


on 


2n+l 2n+l,Q,j 


(35e) 


v (2 ? - X [ 2(n + 1) - 6 ] (2n + 1) a„ ? 2n g„ . (n ) . 

o,i „ L on J 2n+l 2n+l,o,i 

J n=0 ’ 


(35f) 


The superscript (1) indicates the component of the velocity in the 
boundary layer which is in phase with the free stream oscillation, and the 
superscript (2) denotes the out of phase component. 

Figure 1 is a sketch of a typical blunt body and the coordinate system. 
The tangential and normal velocity components are u and v. The Cartesian 
coordinate system (x,y) is also shown in figure 1. The velocity compo- 
nents (u,v) , parallel and perpendicular to the body axis are related to 
(u,v) by 


cos 0 

- v sin 9 , 

(36a) 

sin 9 

+ v cos 9 . 

(36b) 


There is an inherent difficulty in presenting the results of the calcu- 
lations in such a way as to make them intelligible. The reason for this 
difficulty is that the character of the flow, in the boundary layer varia- 
bles, changes drastically with position. In front of the body the flow is 
primarily directed towards the body, so u is very small and v dominates, 
in fact for n large, v « - U 0 . As the flow moves around the body v 
decreases and u increases so that near the midpoint u dominates. For 
this reason it seems most reasonable to present the results in the (x,y) 
coordinate system. In particular, in terms of u, the component of the 
velocity is parallel to the axis of the body. Note, however, that as y is 
varied at fixed x both 5 and n vary. Presentation of the results in 
terms of u also facilitates comparison of the results of these calcula- 
tions with the experimental results of Kachanov, Kozlov, and Levchenko 
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(1978) although, as will be discussed below, the correspondence between the 
conditions of the experiment and of this calculation is not exact. 

Detailed calculations have been carried out for the flow past two types 
of bodies; elliptic and parabolic cylinders. 

3.1 Elliptic Cylinders 

On the surface of an elliptic cylinder the tangential component of the 
potential flow is 

H(5 ) = [l + (b/a)] sn(5;k), (37) 

with a the semi-major axis, b the semi-minor axis, 
sn the Jacobian elliptic function, and 

k? = 1 - (b/a)^ (38) 

The length scale is the semi-major axis, a. The coefficients in the power 

series expansion for H(5), the (a } are given in Table 1, (Cayley, 

zn+ 1 

1895). The radius of convergence of the Blasius series for an elliptic 
cylinder is, apparently, not known. However, the calculations appear to 
converge up to at least the midpoint on the body, that is to 

K = E (k) , (39) 

where E is the complete elliptic integral of the second kind. The 
calculations were ended at this value of 5 because we are primarily 
interested in the flow on the nose. 

Calculations of the velocity field have been carried out for a number 
of bodies with different slenderness ratios (a/b) and frequencies, a. 

Results are given for three bodies; a moderately thick ellipse with a/b = 

5; a thin ellipse, a/b = 10; and a very thin ellipse, a/b = 25, see 
Table 2. For each of these values of a/b calculations were carried out 
for a low frequency a = it / 4 , two intermediate frequencies, a = tr and 
3 tt , and a high frequency, 0=6^. In all cases the Reynolds number was 


14 



taken to be 100.0. The results at other, particularly higher, Reynolds 
numbers appear to have the same form vrtien distances are rescaled to the 
appropriate boundary layer thickness. 


Contours of the zero order steady flow, u , 

0,0 


are shown in Figures 2, 


3, and 4. The format of these figures, and all of the others is the same. 
The body shape is plotted in the lower right hand corner of the figure and 
the flow is from left to right. For each body the semi-major axis is the 
same length. The body shape has not been distorted, however the region 
outside of the body has been stretched in order to show the details of the 
flow in the boundary layer. The coordinate system outside the body is the 
(5 ,n ) system which is shown in Figure 1. The distance 5 = 1 corresponds 
to four of the major units of the scale shown on the edges of Figure 2 and 
all others. The distance H - 1, on the other hand, corresponds to one of 
these major units. Thus, for example, the coordinates of the point in the 
lower left hand corner of Figure 2 are 5=0 and r| = 4, and those of the 
upper right hand corner of Figure 2 are 5 = E(/0. 96) = 1.05050, and n = 
4.2. Note, however, that although the (5 ,0 ) coordinate system is used, 
these are contour plots of u, the velocity component parallel to axis of 
the ellipse. 


The results shown in Figures 2, 3, and 4 show that the zero order in 

e , mean flow has the same general features for bodies of different 

slenderness. Ahead of the body there is a boundary layer merging into the 

decelarating potential flow. Near the midpoint of the body the flow is 

nearly parallel to the body and the boundary layer flow merges into the 

tangential potential flow. In between there is a kind of "bubble" where 

u is small. This extends somewhat in front of and above the forward 

0,0 

• • A 

portion of the body. This region of small u occurs because both u 

o,o 0,0 

and v are small and are at a substantial angle to the x axis. As the 

0,0 

fluid flows around the nose of the body u increases but v decreases 
until near the middle of the body where it is 0(Re“^2). Therefore, on the 
forward portion of the body, but away from the axis, u is increasing with 
5 but is still not large and v is decreasing. Finally, as the body 
becomes more slender, the steady boundary layer tends to become thinner. 
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Figures 5 through 16 are contour plots of the second order steady flow, 

u . Figures 5 to 8 show the results for an elliptic cylinder with 
/ , o 

a/b = 5 for a sequence of frequencies, a = tt/4, tt , 3ir, and 6ir . 

Figures 9 to 12 show similar results for a/b = 10, while the results for 

a/b = 25 are given in Figures 13 to 16. Examination of these figures 

reveals some general patterns in the flow. 

In all cases u. is positive in the lower portion of the boundary 
z , o 


layer and negative in the upper portion. At low frequencies the maximum and 

minimum values of u„ occur at the midpoint of the ellipse. As the 

2,o 

frequency increases the maximum tends to shift towards the front of the body 

and the thickness of the regions of both positive and negative secondary 

mean - flow tend to decrease. However, as the ellipse becomes more slender 

i.e., a/b increases, the region of appreciable u„ tends to be confined 

2,o 

to the near midsection of the body. The absolute magnitude of this secon- 
dary flow also tends to decrease with increasing a. Finally, at low 

frequencies the max (u ) is greater than |min(u )|, but at high fre- 

2,o 2,o 

quencies this is reversed. These results show that the Reynolds stress 
generated by the nonlinear interaction of the fundamental oscillatory flow 
with itself cause the formation of a weak, secondary, shear within the main 
steady boundary layer. At low frequency the maximum shear occurs near the 
middle of the boundary layer. As the frequency increases, the position of 
maximum secondary shear shifts towards the boundary. 

Contour plots of the amplitude and phase of ^ , the fundamental 

oscillation at frequency a, are given in Figures 17 through 28. As above, 

these show results for a/b = 5, 10, and 25 and a = rr/4, ir , 3ir , and 6tt . 

The phase is, of course, arbitrary to within an additive constant. It has 

* 

been assumed that for both the fundamental and the first harmonic, the phase 
is zero at 5=0, n = 4. 

■ 

In all cases the contour maps of the amplitude of u Figures 17a 

o,l, 

to 28a, bear an obvious generic relationship to the corresponding contour 

plots of u for the same values of a/b. Of course, the structure of 

0,0 

the flow does change with the frequency of oscillation. The amplitude of 
the oscillation decreases in front of the body and a region of reduced flow 
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is present above and, generally, in front of the nose of the body. As the 
oscillations move up onto the body, the amplitude grows. In all cases a 
peak in the amplitude distribution appears in the middle to lower part of 
the boundary layer at the midpoint of the ellipse. For a fixed value of 
a/b the location of this maximum of u^ ^ moves towards the boundary as the 

frequency increases, reflecting the development of a Stokes shear layer at 
high frequencies. If the frequency, o, is held fixed, the location of the 
maximum of u^ q also moves deeper into the boundary layer as a/b 

increases, i.e. as the body becomes more slender. This is shown in Table 3, 

in which the location of the maximum of u is tabulated as a function 

o,l 

a and a/b . 

The contour plots of the phase of u Figures 17b to 28b, are all 

o,l 

quite similar. The region of deceleration of the potential flow and the 
development of the boundary layer are obvious, as is the thinning of the 
boundary layer as it develops into a Stokes layer at high frequencies. Just 
back of the nose, the lines of constant phase tend to be parallel to the 
body with a noticeable thickening at low frequency near the midpoint of the 
ellipse. As o increases the region over which the phase shift occurs 
tends to thin. In all cases this region lies deeper in the boundary layer 
than the position of maximum u^ Independently of the value of a/b, 

the total phase shift appears, from the results given in Table 4, to be 
tending to ir/4 as a ■* 00 . 


Finally, Figures 29 through 40 are contour plots of the amplitude and 

/V m w • 

phase of u , the first harmonic of the free stream oscillation. Again 

O, L 

these results show very substantial internal consistency. 


The boundary layer in front of the body is readily apparent. The 
region of nonzero |u | decreases as the frequency increases at constant 
a/b and as the slenderness, a/b, increases at constant frequency. The 
amplitude of u has a maximum in the lower part of the boundary layer. 

O f Jm 

For fixed a/b the position of this maximum moves lower in the boundary layer 
and forward on the body as o increases. At constant frequency the height 
of the maximum remains nearly fixed while moving back toward the midpoint of 
the ellipse as a/b increases. 
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The phase shift of u is generally nearly constant over the region 

O • z 


containing the maximum of |u 


o,2' 


and then goes negative in the outer 


portion of the boundary layer. In this outer region the phase then 
increases towards zero. At low frequencies this outer phase shift occurs 
over a thin layer, but at higher frequencies it is fairly broad. 


The first harmonic of the free stream oscillation is always confined to 

the steady boundary layer and, as o increases, is confined to a thin 

Stokes layer adjacent to the body. As the fluid moves up onto the body, 

away from the nose, the amplitude increases but never becomes large. In 

fact, |u A is never greater than about 5 x 10 - ^. The velocity 
Oj z 

fluctuations due to harmonics of the free stream oscillation thus never 
exceed 5 x 10 -2 g 2 . 


3.2 Parabolic Cylinders 

We next consider the flow past a parabolic cylinder 


y^ = 4 ax, 


(40) 


where a is the distance from the nose to the vertex. Choosing the length 
scale to be a, the tangential component of the potential flow is a power 
series in odd powers of £ = x/a. The coefficients of this series are given 
in Table 5. The solution does not contain any geometric parameters, that 
is, apart from the scaling, is independent of a and is, therefore, a 
similarity solution for all parabolas. 

The radius of curvature at the nose of a parabola is 


r o = 2a. (41) 

Using Van Dyke's (1957) result that the Blasius series diverges for x/r Q 
> 0.62, the solution given here is expected to converge only for £ < 1.24. 

The coefficient aj = 1/4 for a parabolic cylinder, so that 
a = 4 (o)£ /v Q ) . Thus the "effective" frequency in the differential equa- 
tions for the f’s and g's is four times the true frequency of oscilla- 
tion. Of course, a similar scaling applies to elliptic cylinders, but for 
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an ellipse a x = (1 + b/a)~l , which is nearly unity for slender ellipses. 
Because of the factor of four in the value of a for a parabola, 
calculations were carried out for a = m/16, m/4, 3m /4, and 3m /2, so 

that results would be readily comparable to those obtained for the ellipses. 

Contours of u Q Q , the 0(1) mean flow past a parabolic cylinder are 
given in Figure 41. The format of this, and all of the other figures 
showing the results of the calculations for a parabolic cylinder, is similar 
to those figures showing the results for elliptic cylinders. The body is 
shown in the lower right hand corner. Distances normal to the boundary are 
scaled so that one major unit of the scale on the border of the figure 
corresponds to one unit of D, the boundary layer coordinate. Thus the 
point in the lower left hand corner of the figures has the coordinates 
£ = 0.0, r| = 6.0, and the top of the parabola at the right hand side of the 
figures has the coordinates £ = 0.90, n = 0.0. 


A parabola is, in some ways, a rather peculiar blunt body in that it is 
not a closed body, so that there is no trailing edge, and that the thickness 
of the body is a monotonically increasing function of the distance from the 
stagnation point. Therefore, the farther back from the nose, the thicker is 
the region of potential flow which is substantially influenced by the body. 
This effect is apparent in Figure 41. There is a region of reduced flow in 
front of the parabola due to the deceleration of the potential flow and the 
formation of the boundary layer. This region extends far above the axis of 
the body because of the blocking effect of the parabola on the potential 
flow ahead of the body. In the region in front of the nose the contours of 
u q are nearly equally spaced. Just above the parabola, these contours 
tend to line up with the body and the distance between them becomes 
increasingly smaller as the flow moves onto the parabola and the quasi- 
parallel boundary layer flow develops. 


Contours of the steady, second order, streaming flow, U 2 Q > for 
4o = m/4, it, 3tt , and 6m, are shown in Figures 42 through 45. These 
plots show features which are similar to those seen in the contour plots of 



for the flow past an ellipse. Figures 5 through 16. 


In both sets. 
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u_ is positive in the lower portion of the boundary layer and negative in 
2,o 

the upper part. The thickness of the secondary boundary layer decreases 

as the frequency increases and there is a slight tendency for the position 

of the relative maximum and minimum of u^ Q to shift forward towards the 

stagnation point as o increases. However, as o increases the max|u | 

z ,o 

decreases significantly. 

Figures 46 to 49 are contour plots of the amplitude and phase of u , 

o, 1 

the fundamental oscillation in the free stream flow past the parabola. The 

contours of the amplitude of u , show a general similarity to those of 

o, I 

steady flow, u . However, as a increases a Stokes layer develops on 

0,0 

the parabolic cylinder and moves forward towards the stagnation point. Away 

from the stagnation point, the contours of u are similar to those of 

o, 1 

u , for an ellipse if the surfaces are rotated so as to be locally 

parallel. The local maximum of u , seen in the flow past an ellipse is 

o, 1 

not, however, seen for the flow past the parabola. The contour plots of the 

phase of u , also show the development of a Stokes shear layer. The 

0,1 

total magnitudes of the phase shift for u , in Figures 46b through 49b 

o, i 

are quite close to those for the corresponding cases for the elliptic 
cyl inders. 


Contour plots of the amplitude and phase of 


o,2 


the first harmonic 


of the oscillation in the free stream speed, are shown in Figures 50 to 53. 
The results shown in these figures illustrate the development of a second 
order Stokes shear layer within the steady boundary layer as the frequency 
increases. The max|u o increases with £, at all a, at least up to 


the point £ =0.9. It cannot be determined whether or not there is a 
localized maximum in the amplitude, as was found for the elliptic cylinders, 


or whether the amplitude reaches a plateau. Hie problem is that, in order 
to decide this question, the calculations would have to be extended sub- 
stantially beyond £ = 0.9, but, as mentioned above, the series expansion 
probably is divergent for £ > 1.24, and probably converges very slowly for 
£ close to 1.24. 
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DISCUSSION AND CONCLUSIONS 


Perhaps the most general conclusion which can be drawn from the results 
of this study is that in the region of the nose of a symmetric, two dimen- 
sional blunt body at zero angle of attack, the steady plus oscillating flow 
is very similar for a wide class of body shapes. This conclusioon has been 
shown to be true for elliptic cylinders with a/b < 25, and for the parabo- 
lic cylinder. Additional calculations, not reported here, were carried out 
for elliptic cylinders with values of a/b up to 100, with results which 
are very similar to those reported here. In all cases, the flow field in 
the nose region of a two-dimensional blunt body is generic to that of the 
flow in the neighborhood of the stagnation point on a plane wall. 

A general picture of the flow on the forward portion of a blunt body, 

due to a steady plus oscillating free stream, can be sketched. Forward of 

the body there is a region of decelerating potential flow which merges into 

the viscous stagnation point boundary layer. On the stagnation stream line 

the boundary layer thickness is 2 .38^ (Re) - ^^ , with l the scale length 

of the body. If one considers u, the velocity component parallel to the 

body axis, there is a region of small u, above and ahead of the nose, 

where the velocity is reduced and is at a significant angle to the body 

axis. On the body itself, the boundary layer changes slowly with distance 

from the leading edge. The description, so far, applies to u , the 0(1) 

0,0 

steady flow. 

If one next turns to u the fundamental oscillating component of 

o, 1 

the flow, the above description is, with some additions and amendations, 
valid. These are due to the changes with frequency in the oscillating flow. 

At low frequencies, uUl/a.U < — , the flow is essentially quasi-steady; 

° 4 

that is, it is the steady flow scaled by the instantaneous free stream 

speed. As the frequency increases, the oscillating boundary layer develops 

into a Stokes shear layer. This change occurs along, at least, the entire 

forward portion of the body, independent of the position on the body. This 

is not completely obvious when u is examined, but it is obvious if 

o,l 

u q j, the component of the velocity locally parallel to the body is 
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studied. As shown in Appendix B, the development of a Stokes 

layer at high frequency occurs at, and in the neighborhood of, a stagnation 

point. 

This is in marked contrast to the results of Ackerberg and Phillips 
(1972), who concluded that the Stokes layer only develops far downstream of 
the leading edge of a flat plate. This may be due to the fact that there is 
a singularity at the leading edge of the (nonphysical) infinitely thin flat 
plate, or it may be due to the scaling assumptions of Ackerberg and 
Phillips, which, in effect, equate low frequencies with small distance from 
the leading edge and large distance from the leading edge with high frequen- 
cies. In the calculations reported here, the distance from the stagnation 
point and the frequency are independent. 

The oscillations in the boundary layer, which exist at all frequencies 
and at all positions on a blunt body, could be interpreted as an oscillation 
in the steady boundary layer thickness. In the free stream, the mean (time 
averaged) speed is U , and the boundary layer thickness is conventionally 
defined as the height in the boundary layer, <5, where u(6) = 0.99 U q . As 
u increases and decreases at fixed h, because of the oscillation, the 
instantaneous value of 6 would appear to increase and decrease. It is 
quite easy to estimate the amplitude of the apparent oscillation in 6, if 
the velocity oscillations in the free stream are small. Let 

5 = 6 q + «i (t) , (42) 

with S q the time averaged boundary layer thickness, and 5i the oscillat- 
ing component of the boundary layer thickness. By definition, 

<u(6)> = 0.09 <U> = 0.99 U , (43) 

o 

where <> is a time average. Then, if £ << 1, it is easy to see that, to 
lowest order in e , 

cos (to t + 0 ) , (44) 


«1 " 


3u 

e/(_J^_)„ 

9 y y = 6 


0 J 
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where, as above, u is the 

0,0 

parallel to the boundary, and 

y - 6 o- 


0(1) steady boundary layer flow component 

0 is the phase of u , , evaluated at 

o,l 


In the neighborhood of the stagnation point on any blunt, two dimen- 
sional body, 6, scaled by SL , is 

5 o = 2.38 (Re)** 1 7 2 , (45) 

and 


6u 

( — — ) - 2.6 x 10~ 2 . (46) 

/v 

6y y = 5 


So that, in the vicinity of the stagnation point. 


6 l | , scaled by £ , is 


6 ! | “ 38.5 e ( ai Re)- 1/2 


(47) 


and 


l 5 l|/« 


16.2 e( ai ) 


,- 1/2 


(48) 


At low frequencies the flow is very nearly quas i- steady , but at high 

frequencies, there can be some distortion of the velocity profile. If a 

is large, u is essentially constant outside of the Stokes layer which 

o,l 


ome 


has a thickness of 0(6 Inside the layer, u . can cause s 

o J o,l 

distortion of u , while outside the Stokes layer u , only adds a 

0,0 o,l 

constant, in n , time varying increment to u Q Q . 

It may be useful, at this point, to give an idea of what is a low 
frequency and what is a high frequency. From equation (11), 


a ~ ml / (aj U Q ) . 


If, by way of illustration, SL - lm , ai = 1, and Uo = 10 m /sec, then 


23 



a = irf/5 


(49) 


where f is the frequency in Hertz. Thus for 
o < 7r /4 , 

the low frequency region, 
f ~ 5/4 Hertz. 

While, for high frequencies. 


o Z 6ir , f Z 30 Hertz. 


In summary, for this example, frequencies below 1 Hz are low frequencies, 
those greater than 30 Hz are high frequencies, and the range of 1 to 30 Hz 
is the intermediate range of frequencies. 


Considering next the terms of 0(e 2 ), there is both a second order 


steady streaming flow, u. 


and u 


the first harmonic of the free 


Or.* r. n> 

z,o o, Z 

stream oscillation. The steady streaming flow extends throughout the 0(1) 
boundary layer at low frequencies, but at high frequencies, it is largely, 
but not completely, confined to a Stokes layer. However, at high 
frequencies, there is a small portion of u^ q which decays to zero in an 


outer layer, with a thickness of the same order as that of the 0(1) mean 
boundary layer. 


The magnitude of this flow, |u |, increases with £ , and for an 

z , o 

elliptic cylinder has a maximum just in front of the midpoint of the 
ellipse. Such a maximum may, or may not, occur on a parabolic cylinder. 
Because of the limited range of convergence of the series expansion, the 
calculation cannot be reliably extended far enough to decide. 


The first harmonic, u has a structure which is partially similar 

o,2 r J 

to that of the secondary streaming flow, u„ , and partially, to the 

L ,0 
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r 


fundamental, 


o,l' 


At low frequencies it extends throughout the steady 


boundary layer and at high frequencies is confined to a Stokes layer. It 
also develops a maximum amplitude some distance from the stagnation point. 


Both the first harmonic and the secondary streaming flow are present at 
all £ , even at the stagnation point. Away from the stagnation point they 
are somewhat modified, but have the same character as in the neighborhood of 
the stagnation point. In particular, the estimate (see Appendix B) 
that they are bounded by ( e /<*) seems to be valid over an appreciable 
range of 5 , in fact, up to the midpoint on elliptic cylinders and over the 
range of a convergence of the series for flow past a parabolic cylinder. 


The only experiment with which the results of these calculations can, 
apparently, be compared is that of Kachanov, Kozlov, and Levchenko (1978). 
They studied the flow past a flat plate with a nose consisting of two con- 
jugate ellipses with a semi-major axis, a = 50 mm and semi-minor axis, b = 

2 mm on the working side and a = 128 mm and b = 8 mm on the other side. 
Figure 54 shows the cross section of this body from the nose to a distance 
of 20 mm back from the nose. 

In Figure 2 of their paper, Kachanov et al. give contours of (in the 

A. 

notation used here) u in the region from about 10 mm in front of the 

0,0 

nose to about 15 mm behind it, and from the axis of the body to 6 mm 
above. Figure 55 is an enlarged copy, supplied by Dr. Levchenko, of Figure 
2 of their paper. In this figure the contours are labeled with the speed in 
units of meters/second. 

It is obvious that there is a general, qualitative agreement between 
the theoretical results shown in Figures 2, 3, and 4 for the mean flow and 
the experimental results of Kachanov et al. shown in Figure 55. The major 
difference between the theoretical and experimental results is the closed 
contour labeled 6.0 in Figure 55; nothing like this is seen in the results 
of the calculations. The experimental result is rather curious. Either 
there is an absolute maximum within this contour, a "peak", or a local mini- 
mum lies within the contour, so that the contour is the "lip of a volcano." 
In either case the flow speed, as shown in Figure 55, is a maximum outsid e 
of the boundary layer. 
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It is well known that, for a potential flow, the maximum velocity 

occurs on the boundary. The potential flow past an ellipse is a classic 

problem (Milne- Thomson, 1955). Calculations of the potential flow for an 

ellipse with a/b = 25 confirm that the maximum of u occurs on the 

0,0 

boundary at the midpoint of the ellipse. Taking into account the existence 
of the viscous boundary layer, the maximum speed should occur at the top of 
the boundary layer. How is it then possible to account for the existence of 
a maximum in the speed in what is, apparently, a region of potential flow? 

There seems to be three possibilities. The first, and most obvious, is 

experimental error. This does not appear to be likely because it would 

require systematic errors of the order of 10% to 15%. A second possibility 

is that the stagnation point is not on the axis of the body (see Figure 54) . 

If, as seems likely, the stagnation point lies below the axis, then the flow 

must first move through a region of adverse pressure gradient until it gets 

around the nose. Then, for a short distance, the pressure gradient is 

favorable, and then falls rapidly to almost zero just behind the nose. It 

is possible that there is a local flow separation and reattachment, a 

separation bubble, just behind the nose of the upper ellipse. This would 

lift the boundary layer in this region and might account for the closed 

contour of u in the results of Kachanov, Kozlov, and Levchenko. The 

0,0 

third possibility is that the free stream vorticity is not zero. If the 
free stream vorticity is non-zero, perhaps there is a shear in the free 
stream, then, depending on the distribution of free stream vorticity, a 
maximum of the speed could occur virtually anywhere outside of the boundary 
layer. 

There is, of course, a final possibility. These theoretical arguments 
and calculations may be wrong. It may be possible, for some reason, that a 
non- rotational potential flow has a maximum away from the boundaries. Or, 
it may be that there is a viscous effect, not included in these calcula- 
tions, which causes a maximum in the speed outside of the viscous boundary 
layer. 

Kachanov et al. also give experimental results for an oscillatory 
past this body. The oscillatory flow is caused by a vibrating ribbon above 
the axis of the plate and in front of it. As Kachanov et al. indicate. 
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the periodic vortices shed by this ribbon passed above the plate and only 
the tails of the vortices impinged on the nose of the plate. The results of 
the calculations given here are, at best, in rough qualitative agreement 
with the experimental results of Kachanov, Kozlov, and Levchenko. In part, 
their results show an absolute maximum in the amplitude of the oscillatory 
velocity in a small region just above the axis of the plate and just beyond 
the top of the steady boundary layer. They also show a rapid change in the 
oscillatory flow immediately below the region of maximum which was discussed 
above. As was discussed above, either the calculations reported here are in 
error or the experimental results reflect unnoticed effects. In fact, from 
the results given in Figure 1 of Kachanov et al. it appears that the 
oscillation has a nonzero shear in the free stream. In any event, this 
issue requires clarification and would seem to require further work, both on 
the experimental and theoretical aspects of this problem. 

Finally, if the theoretical approach used here and the results are 
correct, it would seem to have value beyond this particular study. The 
general character of the flow, in the nose region of the blunt bodies con- 
sidered here, is very similar to that of the flow in the immediate neighbor- 
hood of the stagnation point. As the flow moves around the nose of the 
body, no essentially new features appear. This suggests that the structure 
of the flow on the forward part of a blunt body due to more general classes 
of disturbances, such as vorticity waves, can be found by studying the flow 
due to these disturbances in the neighborhood of a stagnation point. It 
seems clear, at least for two dimensional flows, that a linearized, in e, 
theory is sufficiently accurate. The results given here show that, for any 
ellipse and for any parabola, the second order terms are bounded by (e/a) 2 
for a £ 1. It also should be noted that the steady flow in the neighbor- 
hood of a general three-dimensional stagnation point is known. This should 
facilitate the theoretical study of disturbances impinging on a body in the 
region of a stagnation point. These approximations, flow near a stagnation 
point, and linearization in e, should be a considerable simplification for 
both analytical and numerical studies of boundary layer receptivity. 
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APPENDIX A 


DIFFERENTIAL EQUATIONS AND BOUNDARY CONDITIONS FOR THE F'S AND G'S 


The function Fj (n , T ) is the solution of the oscillating stagnation 
point flow problem (Grosch & Salwen, Appendix B), The differen- 
tial equations and boundary conditions for the f's and g's are given 
there through terms and will not be repeated here. 


The differential equations and boundary conditions for the f's 


g's given here are those in the expansion of the {F.(ti,t)} for 

J 


J - 3,5, ... ,11 and 


with n < 2. First, we define an operator 


'j ,M \ = 


V 8 J>k ^ dn 


+ f. — - (j + 1) 

,3 l,o,o ^2 l,o,o 


d 

dn 


and 
L by 


+ 


o.J 



(Al) 


where f^ q Q (n) is the 0(1) function in Fi , and is the Hiemenz function, 

and primes denote differentiation with respect to n. Note that in (Al) the 
coefficients are j and j + 1 and that j is the first index of the f 
or g upon which L operates. Next we define a second operator S by 


S(f(n ), g(n);a,S,y) = afg" - Sf'g' + yf"g. 


( A2) 


where f and g are any functions of n, primes, again, denote differen- 
tiation and a, 3, and y are arbitrary constants. 

Then for the terms proportional to we have, at 0(1) 

Lf o = -1 (A3a) 

3,o,o 
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(A3b) 


f, (0) = fl (0) - 0, 
j j O j O J j O j o 


3,0,0 4 


as ri + 


(A3c) 


At 0(e) we have 


Lf, i ' a g'o i = - 2 - S(f . f ; 
3,o,l 3,o,l l,o,l, 3,o,o 


f. ;i,4,3). 


( A4a) 


Lg_ , + of' =— a - S(g., , , f- ; 1, 4,3) , 

6 3,o,l 3,o,l ^ & l,o,l’ 3,o,o’ 


(A4b) 


f 3,o,l (0) f 3,o,l (0) 8 3,0,1 (0) 8 3,0,1 (0) °’ 


( A4c) 


*3,0,! * g 3,o,l " ° 3S n * “ 


(Md) 


And at O(e^) 


if, 9 „ = " 9 - s< f i 9 n 1,4,3) 

3 , 2 , o 2 1 j 2 jO 3,o,o 


(A5 a) 


" 2 f S(f l,o,l’ f 3,o,l’ 1,4,3) + SCg l,o,l ,8 3,o,l ;1,4,3 ^* 


f (0) = f' (0) = 0 
3,2,0 3,2,0 


(A5b) 


f' >0 as n > 
3,2,o 


(A5c) 


and 


Lf - 2ff g' - - - - S(f ,f : 1,4,3) 

3,0,2 3,o,2 2 -I- 3,o,o 


“ 7 £ S(f l,o,l ,f 3,o,l ;1,4,3) “ S(8 1 ,o,1’ 8 3,o, 1’ 1,4,3 ^’ (A6a) 
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Lg, + 20 f ' 

J ) o j 4 3 , 0,2 


- S(g ,f • 51,4,3) 

I , O , Z 3,0,0 


9 ^ S(f l,o,l ,s 3,o,l ;1,4,3) + S(8 l,o,l ,f 3,o,l ;1,4,3) ^ 


(A6b) 


f 3,o,2 C0) * f 3,o,2 (0) = «3,o,2 (0) = 8 3,o,2 (0) = °’ 


9 + „ 9 * 0 as n “• 

3,0, Z 3,0,2 


( A6c) 
(A6d) 


Next the f's and g's proportional to satisfy, first to 0(1) 


L£ 5,o,o » + { C 5> - 8 °5 S(f 3,o,o* l' 1 '^ 

£ 5,o,o (0 > - £ 5,o,o <0) * 0 
f 5 ,o,o + j as 11 * ”• 

Then at 0(e) 

l£ 5 ,o,l “ ° 8 5 ,o,l - 'id + Y C 5 ) ~ S(£ 1 ,o,1* £ 5,o.o: 1 -«>» 

- 4c 5 s££ 3,o,o’ £ 3,o,1 :1 ’ 2,1 ^ 

L8 5,o,1 + o£ 5,o,l “ J a " s(8 l,o,l* £ 5,o,o ;1 " 6 > 5) 

-4C 5 S(£ 3[0j0 ,g 3>0l ; 1,2,1), 

f 5,o,l<°> * f 5,o,l (0> ' 8 5,o,l (0) * 8 5,o,l ( °> ' °> 
f 5 , 0 , 1 " \> 85,o,l * 0 as 11 * “* 

The equations and boundary conditions at O(e^) are: 

“5,2,0 - - { <1 + \ %) - SCfi,2,o, £ 5,o,o;1.6,5) - \ 

* [ s(£ 1,o,1' £ 5,o,1 ;1,6,5) + s<8 1,o,1’ 8 5,o,1 ;1 ’ 6 > 5) ] 

- 8 C 5{ S(£ 3,2,o £ 3,o,o ;1 ‘ 2 - 1) 


(A7 a) 
( A7b) 
(A7c) 


( A8a) 


( A8b) 
( A8c) 
(A8d) 
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— [sC fo n ,.f, . ,i — »1>— ) + S(g q ,,g» 

9 j,o,l 3 ,o,l 2 2 3 ,o,l 3,0,1 


I i ±l) 

2’ ’2 J ; 


f (0) = f' 9 (0) 

3 , 2,0 3 , 2,0 


0; f' 0 as n -»■ «; 

3 , 2,0 


(A9a) 

(A9b,c) 


Lf- 9 - 2ag 

J j O j A 3 j O j 2 


- 2 (1 + 2 °5> S(£ l,o,2 ,£ 5,o,o ;1 ’ 6,5> 


-?I S( 'l,o,l' f i,o,r 1 '‘' i) - S( Sl,o,l’ s 5,o,l ;1 ' 6 ' 5) ] 

- 8C 5( S<f 3,o,o ,£ 3,o,2 :l,2 ’ l) *1 [ S(f 3,o,l' £ 3,o,l ; T’ 1 *? 5 


S(8 3,o,l ,g 3,o,l : 2’ 1 ’2 ) -0 


(AlOa) 


Lg 5,o,2 + 2 ° f 5,o,2 " S(g l,o,2 ,f 5,o,o ;1,6,5) ~ 2 ^ S(f 1 ,o,l ,g 5,0,l ; 

x 1,6,5) + s ^g 1}0jl f 5>O}1 ; 1 » 6 » 5 )] “ 8C 5( SCf 3,o,o ,g 3,o,2 ; 

1,2,1) + lf S(f 3 , 0 ,l jg 3,o, 1 ; T ,1, '2 ) + { S(g 3,o,l ,f 3,o,l ; 


* y,l,j)]}. 


(AlOb) 


f. 9 (0) = f' 9 (0) = g, 9 (0) = g' 9 (0) = 0; 

j j o j Z 3 jO j 2 j y O) 2 j y o ) 2 


^ f' 


5,o,2 6 5 , o , 2 


•♦■0 as n 


(A10c,d) 
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The f's and g' s which are proportional to satisfy the 

equations and boundary conditions, at 0(1) 


Lf_ = 1 (1 + C ) - C S(f ,f ;9 ,24 , 15) , 
7,o,o 7 7 3,o,o 5,o,o 


f_ (0) = f' (0) =0, f' + - as ri + « 
7,o,o 7,o,o 7,o,o g 


At 0(e) the equations and boundary conditions are: 


Lf 7,o,l " ag 7,o,l 


- 2(1 + C ) - S(f . ,f_ ; 1,8,7) 

7 l,o,l 7,o,l 


-c_[s(f ,f ;9 ,24 , 15) + S ( f ,f ; 
7 3,o,o 5,o,l 3,o,l 5,o,o 


* 9,24,15)], 


Lg., , + afl = _ a - s(g, , ; 1 , 8,7) - 

e 7 ,o, 1 7,o,l g 6 l,o,l* 7,o, o’ 


C_[s(f ,g. n . ; 9,24, 15) + s( g ,f ; 

7 3 , 0,0 3 , 0,1 3 , 0,1 3 , 0,0 


* 9,24,15)], 


f, .(0) = f ' ,(0) = g_ .(0) = g! - (0) = 0, 

7,o,l 7,o,l 7,o,l 7,o,l 


fl . + gl . 0 as n +». 

7,o,l 8’ °7,o,l 


While at 0(e 2 ) we find 

Lf 7 ,2 ,o = “ 7 (1 + °5 ) " S(£ l ,2,o ,f 7,o,o ; 1,8,7) " “^^1,0,1* 


f 7,o,l ; 1,8,7) + S(g l,o.l’ 8 7,o,l ;1,8,7) ^ 


(Alla) 

(Allb,c) 

( A12 a) 

( A12b) 

( Al2c) 
(A12d,e) 
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- C 7 {S(f 3,2,o J f 5 J o i o ;9,24,15) + S(f 3,o,o’ f 5,2 s o ;9 > 24 ’ 15) 
+ 4 S(f 3,o,l’ f 5,o,l ;9,24,15) + S(g 3,o,l ,8 5 ,o,r 


x 9,24,15)]}, 


( A13a) 


f 7 _ (0) = V (0) =0, f ' * 0 as n -»■ ». 

7 j 2 jO 7,2,0 7,2,o 


( Al3b , c) 


Lf 7 - 2a g' = - _(1 + C ) - S(f ,f ; 1,8,7) 

7,o,2 7,o,2 2 7 l,o,2 7,o,o 

“ — [ S(f . i»^7 i> >1)8,7) + S (g 1*87 n 1’ 

2 1 » o , 1 7,0,1 l,o,l >,o,l 

X 1,8,7)] - C 7 (S(f 3 _ 0i0 ,£ 5 0i2: 9,24,15) * 
x 9,24,15) + Y [S(f3 j0>1 ,£ 5j0jl ;9,24,15) - s( Sj , 0 j > 85 >Q J 


X 9,24,15)]}, 


( A14a) 


Lg_ . + 20 fl , 
& 7,o,2 7,o,l 


= - s(g f ; 1 , 8,7) - ^ [s(f , ,g 7 , ; 

X, 0,2 7, 0,0 2 1,0,1 / , 0 j 1 


X 1,8,7] + S(g f , ; 1,8,7) - C {s (f ,g ; 

1,0,1 7,0,1 7 j,o,o 5,o,2 

x 9,24,15) + S(g f ;9 ,24, 15) + i[s(f - ,g, o . ; 

J j O j 2 3)0)0 2 J , O , 1 3)0)1 


9,24,15) + S(g f ;9,24,15)]>, 

J ,o , 1 5,0,1 


( A14b) 


f 7 n ,(0) = f ' ,(0) = g 7 _(0) = g' (0) = 0, ( A14c) 

7,0,2 / , O ,2 / , 0,2 7,0,2 


f 7 „ 9 * «7 o 9 * 0 38 n * “■ 

7,0,2 7,o,2 


( A14d) 
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For the terms proportional to the 0(1) equations and boundary 

conditions are: 


Lf 


9,o,o 


(1 + C Q - + 7r Cq « ) "T7T C- 1 , f-j | 

9,1 2 9,2 10 9,1 3,o,o’ 7,o,o* 


96,320,224) - 18c 9> 2 S( f 5 ,o ,o 5 f 5,o,o ; 2* 2 J ’ 


( A15 a) 


f. (0) = fl (0) = 0, fl ♦ 1/10 as h -*• «. 
9,o,o 9,o,o 9,o,o 


(A15b ,c) 


The 0(e) equations and boundary conditions are: 


Lf - Og' = -2(1 + C n + C n n ) - S(f, , ,f ; 1,10,9) 

9,o,l e 9,o,l 9,1 9,2 l,o,l 9,o,o 

- C Q Js(f, ,f_ . ;96,320,224) + S(f, _ . , f 7 ; 

10 9,1 L 3,0,0* 7,o,l’ ’ 3,0,1 7,o,o 

x 96,320,224)] - 18c S(f ,f :1,2,1) (A16a) 

9,2 5 ,o ,o 5,o,l 

Lg + Of* , “ tk a ~ S(g, , ,f Q ; 1, 10, 9) - ytp C Q . 

9,o,l 9,o,l 10 °l,o,l ’ 9,o,o’ 10 9,1 

x [s(f g ;96 ,320 ,224) + S(g f ; 96 ,320 ,224) ] 

3 j O ^ O / j O j 1 j j O j 1 / j O j o 


- 18C S(f g ; 1,2, 1) , 

9 j ^ I) jO j o !) jO j 1 


(Al6b) 


f q ,(0) = fl .(0) = g Q .(0) = g! .(0) = 0, 

9,o,l 9,o,l 9,o,l 9,o,l 


( Al6c) 


f 9,o,l * "nr* g 9,o,l 


-*■ 0 as n <*> . 


(A16d, e) 


The 0(e z ) equations and boundary conditions are: 


b f Q _ = - 1(1 + C Q . + - C ) - S(f. . f z 

9,2,0 2 9*1 2 9,2 l,2,o 9,o,o 
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1,10,9)] [s ( fj, o>1 .*9, o ,l :l - 10 * 9) + S(g l.o.l’ g ».«. li 


* 1.10,9)] -io S,lf«3. l . O *7. O . O ‘ ,6 ' Sa0 * m> + «W 
x f 7 2 ,o’ , 6 ’ 320 ’ 224) *^ S<f 3,o,l’ £ 7,o,l ;96 ’ 320 ’ 224) 

* S( 83,o,l’ g 7,o,l i96>320 ’ 224)] J ' 18C 9,a( S(, 5,2,o> f 5,o,o ; 

x 1,2,1) ♦ it S(£ 5,o,l f 5,o,l ; 2’ 1 '2‘ ) * S(s 5,o,r g 5,o,l’ 


, 1 , 1 )]} 


9,2,0 


2**2 


(0) - f; >2>o (0) - 0, f; >2>0 ^ 0 as D ♦ 


(Al7a) 

(Al7b,c) 


L£ ( 


l>0i2 - ^9,0,2 * * T (l + °9,1 * I C 9,2> - S£f l,o,2* f 9,o,o’ 1,10,9) 


-T [ S(£ 1,o,1- £ 9,o,1 ;1 > 10 - 9> ' S(g l,o,l’ 8 9,o,l 


;1,10,9) ] 


Lg< 


-T0 C 9,l{ S ' £ 3,o, O £ 7, O ,2 ; 94 - 320 > 224)+S(£ 3,o,l' 

x f ;96,320,224) + |[ S(£ 3 ,o,l ’ £ 7 ,o, 1 ;96, 32 °’ 224> 

/ ,o , o 

- S(g 3 Q ! >87 j0j l’ 96,320 ’ 224 ^ 1 ” 18C 9,2{ S(f 5,o,o’ 5,o,2 
x 1,2,1) + |[ s(f 5,o,l £ 5,o,l ; T 1 ^ ' S(8 5,o,l ,g 5,o,r 

x ■? ,1, T > ^}’ 

,, 0t2 * ^,o,2 • - S£g l,°,2’ f 9,o,o’ 1,10 ’ 9) "I t S<f l,o,l’ g 9,o,l ; 


(A18a) 


36 



X 1,10,9) - S<8 ljOjl .f 9 , o>1 >l,10.9)] - To C 9,l S(f 3,o,o- 

* g 7,o,2> 96 ’ 320 ’ 224) * S(g 3jOj2 ,f 7jOjo! 96,320,224) + I 

* t S(f 3,o,l' g 7,o,l’ 96,320 ’ 224) + S(s 3,o,l’ f 7,o,l ;96 ’ 320, 

X 224)]} - 18C 9 i 2( S<f 5,o,o’ g 5,o,2 ; l ,2 ’l> + S<f 5,o,l’ s 5,o,l ; 


X 1 


, 2 , 1 )}, 


(Al8b) 


f 9,o,2 (0) f 9,o,2 (0) g 9,0,2 (0) g 9,o,2 (0) “ °» 


( A18c ) 


f 9,o,2 * g 9,o,2 - ° 38 n ~ 


(A18d) 


Finally, for the terms proportional to 5^ the differential equations 
and boundary conditions are, at 0(1): 


Lf.. - - (1 + C.. , + C 1t „) - C.. ,S(f- ,f Q ; 10,40,30) 
ll,o,o 11,1 11,2 11,1 3,o,o 9,o, o’ 


- C.. -S(f- ,f_, ; 20,48, 28) 

11,2 5,o,o 7,o,o 


( A19a) 


f,, (0) = f' (0) = 0, 

11 , 0,0 11 , 0,0 ’ 


(A19b) 


f \ , -*■ as r| -»■ 

ll,o,o 12 


(A19c) 


The 0(e) differential equations and boundary conditions are: 


Lf.. . - og' . = -2(1 + C.. , + C.. 0 ) - S(f. ..f.. ; 

ll , o,l ll , o,l 11,1 11,2 l , o,l ll , o,o 


1,12,11) - C llil [ S (f 3jO _ <) ,f 9iOjl ;10,40,30) 


S(f, O i * f o o i j 10,40, 30) ] - C [S(f f ; 

j j o 1 1 7 jO y 1 11)2 JjO jO 7)0)1 
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And, 


x 20,48,28)] + S(f ,f ; 20 ,48,28) ] , (A20a) 

3)0)1 / )0 • O 


Lg i l o l + af i'i n 1 = TT a “ S (g. 1 ,f. . ; 1,12,11) 

-L i. j O j .L -L i. ^ O ) -L 1 4 1)0)1 11 ) O ) o 


-C.. ,[s(f. ,g_ . ; 10 ,40,30) + S(g ,f ; 

11,1 3,o,o ° 9 , o,l 3,o,l 9,o,o 

x 10,40,30)] - C [ S ( f ,g_ ;20 ,48 ,28) 

11,2 o,o,o /,o,l 

S(8 5,o,l ,f 7,o,o ;20,48 ’ 28) ]’ 


( A20b) 


f.. (0) = f ' , (0) = g (0) = g' (0) = 0, 

11 , 0,1 11 , 0,1 11 , 0,1 11 , 0,1 


( A20c) 


f' 


ll,o,l T2’ B ll,o,l 


■* 0 as n -► 00 . 


(A20d,e) 


9 

lastly, the 0 (e^ ) equations and boundary conditions are: 


Lf n,2,0 = - 2 (l + C ll,l + C ll,2 ) - S(£ l,o,2 ,£ ll,o,o ;l,12,11) 

“ 2 f S(£ l,o,l' f ll,o,l !l,12,11) * S<8 l,o,l’ g ll,o,l ; 

X 1,12,11)] - C Utl {8(f J>2tO .f 9iOiO il0.*0.30) ♦ S(f . 


£ 9,2,0' l0 ' 40 ’ 30) ^[ S(f 3,o,l' f 9,o,l’ 10 ' 40 ' 


30) 


+ S(g 3,o,l- g 9,o,l- 10 - 40 - 30) ]}- C 11,2{ S(£ 5,2 ,o’ £ 7,o,o ; 
x 20,48,28) * S(f 5iOo ,f 7>2>o ;20,48,28) + £ [s(f 5-0il , 

X f ;20,48,28) + S(g ,g ; 20 ,48 ,28) ] } , (A21a) 

7 , 0,1 0 , 0 , 17 , 0,1 ’ 


f 11 0 n (0) = f il 9 _ (0) = 

11 , 2,0 11 , 2,0 


(A2 lb) 
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11,2,0 * 0 as n ♦ «. 


(A21c) 


Lf 


2a 8i lj0i2 * (1 ♦ C u>1 t C u>2 ) - S<f 1>o 2 ,f llj0 o ; 

xi,i2,n) -lf s <f lj0>1 .f llj0il! i > i 2 ,iu - s(g 

1,0,1’ 

*(* 3>OiO .* 9>Ot2i 10. W .30) 

+ S(£ 3.o. 2 . f 9 , o , o U0,40,30) *». ^^,^^0.40.30) 
S(S 3>o,l >g 9,o,l ;I0 > 40 > 3 °)] - C 71 0 s(f c f f 


1J »2 5,o,o’ f 7,o,2’ 

x 20,48,28) + SCf 5jO)2 ,f 7jO>o . 20j48}28) + l [g(ft 
x f. 


Lg 


2 L ”'^,o,l 

7,o, 1’ 20 ’ 48 »28) - S( S 5jOjl ,g 7jO ^ ;20j48j28) j ? (A22a) 

U,o,2 + 2af ll, 0 ,2 = " S(g l,o,2’ f H, 0 , o ;l,12,ll) 


30) 


2 ^ s(f i,o,i >8 n >0> i;M2,ii) + s(g 

1 i o , 1 

U, °’ l,I,12 ’ U) ] * C H,l( S(t 3, < , )O .8 9>o 2 ;l0 > 40, 


+ S(8 3,o, 2 .'9, o , o U0,40,30) 4t«* s . o , 1 .« 9>Ofl ,10.*0.J0) 
* S<8 3,o,l* f 9,°,I :l0 > 40 .30)]} - c n <2 { s Cf 5j o i0 . g7j0j2 ; 

« 20.«.«> * « w !20 .48, 2 8) + j [s(£ 

c J»o,l 

X8 7>o , l!2 0,48, 2 8) + s( g5iO>l , f7jOiii20>48i28)])j aj2b) 
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( A22c) 


£ 11,o, 2 (0) ' £ ll,o,2 (0) - Sll,o,2 <0) - 8 [l,o,2 (0) " °* 


£ ! i « 9 + g! , „ „ o as n -> ». 

11,0,2 ll,o,2 


40 



REFERENCES 


Ackerberg, R.C. & Phillips, J.H. 1972, J. Fluid Mech. 51 , 137. 

Cayley, A. 1895, "Elliptic Functions" Second Ed. Reprinted Dover Pub., New 
York, 1955. 

Carrier, G.F. & DiPrima, R.C. 1957, J. Math. & Phys. ^5^, 359. 

Gibson, W.E. 1957, Ph.D. Math. Thesis, Mass. Inst. Tech. 

Glauert, M.B. 1956, J. Fluid Mech. _1, 97. 

Grosch, C.E. & Salwen, H. 1980, Proc. 15th Int. Cong. Theor. Appl. Mech., 
155, Spr inger-Ver lag, Berlin. 

Harvey, W,D. & Bobbitt, P.J. 1981, AIAA-81-1225. 

Kachanov, Yu. S. , Kozlov, V.V. & Levchenko, V. Ya. 1978, Isv. Akad. Nauk 
SSSR, Mekh. Zhid. Gaza 5, 85. 

Klebanoff, P.S., Tidstrom, K.D., & Sargent, L.M. 1962, J. Fluid Mech. J^2, 1. 
Lam, S.H. & Rott, N. 1960, Cornell Univ. GSAE Rep. AFOSR TN-60-1100. 
Lighthill, M.J. 1954, Proc. Roy. Soc. , A, 224 , 1. 

Lighthil 1, M.J. 1978, J. Sound & Vib. 6_1, 391. 

Lin, C.C. 1956, Proc. 9th Int. Cong. Theor. Appl. Mech., 136, Springer- 
Verlag, Berlin. 

Mack, L.M. 1977, JPL Pub. 77-15, Jet Prop. Lab., Cal. Inst. Tech. 

Milne- Thomson, L.M. 1955, "Theoretical Hydrodynamics", 3rd Ed. MacMillan 
Comp., New York. 

Moore, F.K. 1951, NACA TN 2471. 

Morkovin, M.V. 1978, AGARDograph No. 236, NATO, Paris. 

Obremski, H. J., Morokovin, M.V., & Landahl, M. 1969, AGARDogroph No. 134, 
NATO, Paris. 

Riley, N. 1975, SIAM Rev. 17, 274. 

Rott, N. 1956, Quart. Appl. Math. _13, 444. 

Rott, N. & Rosenzweig, M.L. 1960, J. Aerospace Sci. 27, 74. 


41 



Salwen, H. , Kelly, K.A. & Grosch, C.E. 1980, Bull. Am. Phys. Soc. 25_ t 1085. 
Salwen, H. & Grosch, C.E. 1981, J. Fluid Mech. 104 , 445. 

Sarma, G.N. 1964, Proc. Camb. Phil. Soc. 60 , 137. 

Schlichting, H. 1979, "Boundary Layer Theory" McGraw-Hill Comp. New York. 

Smith, A.M.O. 1956, Proc. 9th Inti. Cong. Theor. Appl. Mech. Springer- 
Verlag, Berlin. 

Stuart, J.T. 1966, J. Fluid Mech. 24 , 673. 

Van Dyke, M. 1956, Proc. 9th Inti. Cong. Theor. Appl. Mech. 318, Springer- 
Verlag, Berlin. 

Van Ingen, J.L. 1956, Rept. UTH-74, Dept. Aero. Eng. Univ. of Tech. Delft. 
Watson, J. 1959, Quart. J. Mech. Appl. Math. 12, 175. 


42 



Table 1. Coefficients in the power' series expansion of the potential velo- 
city function, H(? ) , for flow past an elliptic cylinder. 


2n + 1 


1 

3 

5 

7 

9 

11 




1 


-(1 + k 2 ) / 3 ! 

(1 + 14k 2 + k 4 )/5! 

- (1 + 135k 2 + 135k 4 + k 6 ) / 7 ! 

(1 + 1228k 2 + 5478k 4 + 1228k 6 + k 8 )/9! 

(1 + 11069k 2 + 165826k 4 + 165826k 6 + 11069k 8 + k 1 0 )/ 1 1 ! 


Table 2. Parameters of the elliptic cylinders for which results are given. 


a/b 

k 2 

E(k) 

5 

0.96 

1.05050 

10 

0.99 

1.01599 

25 

0.9984 

1.00329 
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Table 3. The 

position in 

the boundary 

layer of the 

maximum of u, 

1,0 

a 

a/b = 

5.0 

10.0 

25.0 

tt /4 


2.2 

2.1 

2.0 

IT 


1.9 

1.8 

1.7 

3ir 


1.2 

1.1 

1.05 

6ir 


0.9 

0.8 

0.75 


Table 4. 

Total 

phase shift. 

in degrees, of 

u 

1,0 

a 

a/b = 

5.0 

10.0 

25.0 

ir /4 


11.2 

11.9 

12.0 

ir 


30.0 

30.0 

30.0 

3tt 


42.0 

40.0 

40.0 

6ir 


42.0 

42.0 

40.0 
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Table. 5. Coefficients in the power series expansion of the potential 
velocity function, H(5), for flow past a parabolic 
cyl inder. 


2n + 1 


1 

3 

5 

7 

9 

11 


1/4 

- 1/2 2 .3! 
22/2 4 .5! 
-1168/2 6 .7! 
113536/2 8 .9! 

- 1752 1024/2 1 0 . 11 ! 
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Figure 1. 


Coordinate system. 



Figure 2. Contours of the zero order steady flow, u Q Q . Contours from 0.0 

to 0. 90 with an interval of 0. 05. The body is an ellipse with 
a/b « 5. 


UlLLIl 




Figure 3. Contours of the zero order steady flow, u Q 0 * Contours from 0.0 

to 0. 80 with an interval of 0. 05. The body is an ellipse with 
a/b - 10. 



■Ill 



Figure 4. Contours of the zero order steady flow, u_ rt . Contours from 0.0 

O j U 

to 0.74 with an interval of 0.04. The body is an ellipse with 
a/b « 25. 




Figure 5. Contours of the second order steady flow, Q . Contours from 

-1. 6 x 10“ 2 to 5. 2 x 10“ 2 with an interval of 4. 0 x 10” 3 . The 
labels are scaled by 10\ The body is* an ellipse with a/b « 5 
and a m tt/4. 


I 




Contours of the second order steady flow, u« • Contours from 

-2. 1 x 10“ 2 to 2. 7 x 10 2 with an interval of 3.0 x 10-3. The 
labels are scaled by 10 4 . The body is an ellipse with a/b ■ 5 
and a » ir. 


I 














Figure 9. Contours of the second order steady flow, U 2 Q . Contours from 

-1.2 x 10 -2 to 4.5 x 10” 2 with an interval of 3.0 x 10” 3 . The 
labels are scaled by 10^. The body is an ellipse with a/b ■ 10 
and 0 ■ x/4. 
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Figure 10. Contours of the second order steady flow, u 2t0 » Contours from 

-2. O x 10 -2 to 2.4 x 10” 2 with an interval of 2.0 x 10~ 3 . The 
labels are scaled by 10 4 . The body is an ellipse with a/b - 10 



Figure 11. Contours of the second order steady flow, U 2 0 » Contours from 

-1.4 x 10~ 2 to 1.0 x 10“ 2 with an interval of 1.0 x 10”^. The 
labels are scaled by 10 1 *. The body is an ellipse with a/b ■ 10 




Figure 12. Contours of the second order steady flow, “ 2 , 0 * Contour8 from 

-B.0 x 10" 3 to 4. 8 x 10*~ 3 with an interval of 8.0 x 10 4 . The 
labels are scaled by 10 5 . The body is an ellipse with a/b - 10 
and o ■ 6ir . 
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Figure 13. Contours of the second order steady flow, u~ . Contours from 

i j o 

-1.2 x 10~ 2 to 4.2 x 10- 2 with an interval of 3.0 x 10” 3 . The 
labels are scaled by 10^ . The body is an ellipse with a/b ■ 25 
and a * tt/4. 



Figure 14. Contours of the second order steady flow, Contours ^ rom 

-1.8 x 10“ z to 2. 2 x 10” 2 with an interval of 2.0 x 10“ 2 . The 
labels are scaled by 10 4 . The body is an ellipse with a/b * 25 





Figure 15. Contours of the second order steady flow, U 2 Contours from 

-1.3 x 10" 2 to 8.0 x 10“3 with an interval of 1.0 x 10”^. The 
labels are scaled by 10 1 *. The body is an ellipse with a/b ■ 25 
and o * 3if. 


A 

Figure 16. Contours of the second order steady flow, u£ Q . Contours from 

-7.7 x 10” 3 to 4.2 x 10” 3 with an interval of 7.0 x 10” 1 *. The 
labels are scaled by 10^. The body is an ellipse with a/b ■ 25 
and o ■ 6tt. 










Figure 18a. Contours of the amplitude of u_ -i , the fundamental oscilla- 

O j J- 

tion. Contours from 0.0 to 0.96 with an interval of 0.06. The 
body is an ellipse with a/b ■ 5, and a m ir. 





Figure 18b. Contours of the phase, in degrees, of u Q <j_, the fundamental 

oscillation. Contours from -28. 0° to 2. 0° with an interval of 
2.0°. The body is an ellipse with a/b « 5, and a m 5. 







Figure 19b. 


Contours of the phase, in degrees, 

oscillation. Contours from -38. 0* 
2.0°. The body is an ellipse with 


of u rt i , the fundamental 

O, J. 

to 4. 0° with an interval of 
a/b - 5, and a ” 3n. 









Figure 21a. Contours of the amplitude of u Q j, the fundamental oscilla- 
tion. Contours from 0.0 to Ov 85 -with an interval of 0.05. The 
body is an ellipse with a/b ■ 10, and a * tt/4. 



Figure 21b. Contours of the phase, in degrees, of u the fundamental 

0,1 

oscillation. Contours from -11.2° to 0.7° with an interval of 
0.7°. The body is an ellipse with a/b - 10, and a m w/4. 








Figure 22b . Contours of the phase, in degrees, of u_ •. , the fundamental 

oscillation. Contours from -28. 0° to 2. 0° with an interval of 
2.0°. The body is an ellipse with a/b = 10, and o ■ it. 





Figure 23a. 


Contours of the amplitude of . , the fundamental oscilla- 

tion. Contours from 0.0 to 0.90 with an interval of 0.05. The 
body is an ellipse with a/b ■ 10, and a ■ 3tt. 



tn 



Figure 23b. Contours of the phase, in degrees, u Q 1 , the fundamental 

oscillation. Contours from ~36. 0° to 4. 0° with an interval of 
2. 0°. The body Is an ellipse with a/b « 10, and o m 3tt. 



Ov 


Figure 24a. Contours of the amplitude of vL -i , the fundamental oscilla- 

Oj J. 

tion. Contours from 0.0 to 0.90 with an interval of 0.05. The 
body is an ellipse with a/b * 10, and o * 6tt. 





Figure 24b. Contours of the phase, in degrees, of u the fundamental 

oscillation. Contours from -40.0° to 2.0° with an interval of 
2.0°. The body is an ellipse with a/b ** 10, and a * 6tt. 




Figure 25a. Contours of the amplitude of u Q the fundamental oscilla- 
tion. Contours from 0.0 to 0.80 with an interval of 0.05. The 
body is an ellipse with a/b * 25, and o * r/4. 



Figure 25b. Contours of the phase, in degrees, of u_ the fundamental 

o, J. 

oscillation. Contours from -12. 0° to 0. 0° with an interval of 
0.8°. The body is an ellipse with a/b ■ 25, and a ■ ir/4. 




Figure 26a. Contours of the amplitude of u Q the fundamental oscilla- 
tion. Contours from 0.0 to 0.85 with an interval of 0.05. The 
body is an ellipse with a/b ■ 25, and o m ir. 


Figure 26b. Contours of the phase, in degrees, of u the fundamental 

O, X 

oscillation. Contours from -28. 0° to 2. 0° with an interval of 
2.0°. The body is an ellipse with a/b ■ 25, and a * it. 







Figure 28a. Contours of the amplitude of u rt ■* , the fundamental oscilla— 

o, i 

tion. Contours from 0.0 to 0.80 with an interval of 0.05. The 
body is an ellipse with a/b ■ 25, and a m 6ir. 





Figure 29a. Contours of the amplitude of u , the first harmonic of the 

r , ^ 

free stream oscillation. Contours from 0.0 to 4.4 x 10" ^ with 
an interval of 2.0 x 10~ 3 • The labels are scaled by 10** . The 
body is an ellipse with a/b - 5, and a = it /4 . 





Figure 29b. Contours of the phase, in degrees, of u rt <>, the first 

harmonic of the free stream oscillation. Contours from 
~200.0° to 140.0° with an interval of 20.0°. The body is an 
ellipse with a/b « 5, and a = ir/4. 


Figure 30a. Contours of the amplitude of u Q 2» ff rst harmonic of the 

J r> 

free stream oscillation. Contours from 0.0 to 2.0 x 10“ * with 
an interval of 1.0 x 10“ 3 . The labels are scaled by 10 4 . The 
body is an ellipse with a/b ■ 5, and o ■ r. 
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Figure 30b. Contours of the phase, in degrees, of v the first 

harmonic of the free stream oscillation. Contours from 
-240.0 to 100.0° with an interval of 20.0°. The body is an 
ellipse with a/b = 5, and a » ir . 




Figure 31a. Contours of the amplitude of % } 2 > the f *- rst harmonic of the 

free stream oscillation. Contours from 0.0 to 4.6 x 10 with 
an interval of 2.0 x 10 -4 • The labels are scaled by 10 . The 
body is an ellipse with a/b = 5> and o — 6 tt . 
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Figure 31b. Contours of the phase, in degrees, of u 0j2 > 

harmonic of the free stream oscillation. Contours from 
-260.0° to 80.0° with an interval of 20.0°. The body is an 
ellipse with a/b ■ 5, and o * 3n. 



Contours of the amplitude of u Q 2 » the first harmonic of the 

free stream oscillation. Contours from 0.0 to 4.6 x 10” ^ with 
and interval of 2.0 x 10" 1 *. The labels are scaled by 10 5 . Th 
body is an ellipse with a/b *= 5, and o = 6n. 





Figure 33a. Contours of the amplitude of u Q ^ ^ tlie first harmonic of the 

free stream oscillation. Contours from 0.0 to 4.0 x 10 2 with 
an interval of 2.0 x 10" 3 . The labels are scaled by 10 . The 
body is an ellipse with a/b - 10, and o - n/4. 





Figure 33b. Contours of the phase, in degrees, of u 0>2 * the first 

harmonic of the free stream oscillation. Contours from 
-200.0° to 140.0° with an interval of 20.0°. The body is an 
ellipse with a/b » 10, and a * tt/4. 




Figure 34a. Contours of the amplitude of u Q 2 * the first harmonic of the 

free stream oscillation. Contours from 0. 0 to 1. 8 x 10 ^ 
with an interval of 1.0 x 10“ 3 . The labels are scaled by 10 4 . 
The body is an ellipse with a/b a 10, and o * it. 
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Figure 35a. Contours of the amplitude of u 0 2* fc he harmonic of the 

free stream oscillation. Contours from 0. 0 to 7. 2 x 10~ 3 with 
an interval of 4.0 x 10” 4 . The labels are scaled by 10^. The 
body is an ellipse with a/b * 10, and a ■ 3ir. 
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Figure 36a. Contours of the amplitude of u q 2> the first harmonic of the 

free stream oscillation. Contours from 0.0 to 4.0 x 10“ 3 with 
an interval of 2.0 * 10" 4 . The labels are scaled by 10 . The 
body is an ellipse with a/b = 10, and 0=6^. 





Figure 36b. Contours of the phase, in degrees, of u , the first harmonic 

£ ^ 
of the free stream oscillation. Contours from -280.0° to 60.0° 

with an interval of 20.0°. The body is an ellipse with a/b = 10 
and a = 6^ , ’ 
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Figure 37a. Contours of the amplitude of u o,2* ^ r8t: harmonic of the 

free stream oscillation. Contours from 0.0 to 3.6 x 10“ 2 with 
an interval of 2.0 x 10“ 3 . The labels are scaled by 10 4 . The 
body is an ellipse with a/b » 25, and a - n/4. 
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Figure 38a. Contours of the amplitude of u Q 2 > the first harmonic of the 

free stream oscillation. Contours from 0. 0 to 1. 6 x 10“ 2 with 
an interval of 1.0 x 10~ 3 . The labels are scaled by 10 4 . The 
body is an ellipse with a/b ■ 25, and o •= r. 



A , 

Figure 38b« Contours of the phase, in degrees, of u the first harmonic 

o , z 

of the free stream oscillation. Contours from -”240.0° to 100.0 
with an interval of 20.0°. The body is an ellipse with a/b =* 25, 
and o =.ir . 
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Figure 39a. Contours of the amplitude of u Q 2 > the first harmonic of the 

free stream oscillation. Contours from 0. 0 to 6. 4 x 10“ 3 with 
an interval of 4. 0‘ x 10~ 4 . The labels are scaled by 10 5 . The 
body is an ellipse with a/b ■ 25, and a - 3tt. 
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Figure 39b. Contours of the phase of u Q 2» the first harmonic of the free 

stream oscillation. Contours from ~260. 0° to 60. 0° with an 
interval of 20.0°. The body is an ellipse with a/b • 25, and 
a m 3ir. 
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Figure 40a. Contours of the amplitude of u Q 2 * first harmonic of the 

free stream oscillation. Contours from 0.0 to 3.4 x 10“ ^ with 
an interval of 2.0 ,x 10 _lf . The labels are scaled by 10^. The 
body is an ellipse with a/b “ 25, and a m 6t. 



Figure 40b. Contours of the phase, in degrees, of G 9 » the first harmonic 

Oj z 

of 'the free stream oscillation. Contours from -280.0* to 60.0 
with an interval of 20.0°. The body is an ellipse with a/b " 25, 
and o = . 


the fundamental oscilla- 


Figure 41. 


Contours of the amplitude of u . , 

0,1 

t ion of the flow past the parabola with a = ir/16. 
from 0.0 to 0.22 with an interval of 0.01. 


Contours 
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Figure 42. Contours of the second order steady flow, Q , past the 

parabola with a = tt / 16 . Contours from -1.4 x 10 -3 to 1.05.x 
10- 2 with an interval of 7.0 x 10*^. The labels are scaled by 
10 4 
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Figure 43. Contours of the second order steady flow, U 2 0 , past the 

parabola with a * tf/ 4. Contours from -3.6 x 10~ 3 to 6.6 x 10~ 3 
with an interval of 6.0 x lO" 4 . The labels are scaled by 10 5 . 




Figure 44. Contours of the second order steady flow, U 2 Q , past the 

parabola with a = 3n/4. Contours from -3.0 x 10” 3 to 2.4 x 
10” 3 with an interval of 3.0 x 10 _l+ . The labels are scaled by 
10 5 . 
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Figure 45. Contours of the second order steady flow, U 2 0 » past the 

parabola with a *= 3ir/2. Contours from -2.0 x 10“ 3 to 1.2 x 
10” 3 with an interval of 2.0 x 10" 4 . The labels are scaled by 
10 s . 


I 
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Figure 46b. 


Contours of the phase, in degrees, of u^, the fundamental 
oscillation of the flow past the parabola with o - tt/16. Contours 
from - 8 . 5 ° to 0.0° with an interval of 0.5°. 
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oscillation of the flow past the parabola with o =* ir/4. Contours 
from -24.0° to 1.0° with an interval of 1.0°. 
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Figure 48a. Contours of the amplitude of u rt , , the fundamental oscilla- 

0,1 

tion of the flow past the parabola with a ■ 3u/4. Contours 
from 0.0 to 0.22 with an interval of 0.01. 
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oscillation of the flow past the parabola with o = 3tt/ 4. Contours 
from -34.0° to 2.0° with an interval of 2.0° 




Figure 49a. Contours of the amplitude of %,1* ^ un ^ ainent:a ^ oscilla- 

tion of the flow past the parabola with o m 3 tt/ 2. Contours 
from 0.0 to 0.22 with an interval of 0.01. 




Figure 49b. 


Contours of 

oscillation 
from -38.0° 


# A, 

the phase, in degrees, of u 


o,l* 


the 


of the flow past the parabola with o 
to 2.0° with an interval of 2.0°. 


fund ament a 1 
* 3 it/ 2. Contours 
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Figure 50a. Contours of the amplitude of u Q 2 » th e fi rB ^ harmonic of the 

free stream oscillation of the flow past the parabola with 
o * tt/ 16. Contours from 0.0 to 9.6 x 10“ 3 with an interval of 
6.0 x 10 " 1+ . The labels are scaled by 10 5 . 
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Figure 50b. 


Contours of the phase, in degrees, of u 9 , of the first harmonic 

o , z 

of the free stream oscillation of the flow past the parabola with 
a - tt/ 16, contours from -260.0° to 140.0° with an interval of 


20 . 0 ° . 




Figure 51a. Contours of the amplitude of u 2 , the first harmonic of the 

o, 

free stream oscillation of the flow past the parabola with 
a - ir/4. Contours from 0. 0 to 4. 6 x lCT 3 with an interval of 
2. 0 x 10 -1 *. The labels are scaled by 10 3 .* 
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Figure 51b. Contours of the phase, in degrees, of u q 2 , of the first harmonic 

of the free stream oscillation of the flow past the parabola 
with o = tt/ 4, contours from to -240.0“ to 100.0* with an interval 
of 20.0*. 


rnmrrr 




Figure 52a. Contours of the amplitude of u Q first harmonic of the 

free stream oscillation of the flow past the parabola with 
a m 3tt/ 4. Contours from 0. 0 to 1. 8 * 10 ^ with an interval of 
1.0 x 10"^. The labels are scaled by 10 5 . 
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Figure 52b. Contours of the phase, in degrees, of u , of the first harmonic 

of the free stream oscillation of the flow past the parabola 
with a = 3tt/4, contours from -200.0° to 80.0° with an interval of 
20.0°. 




Figure 53a. Contours of the amplitude of u q ^ the first harmonic of the 

free stream oscillation of the flow past the parabola with 
a = 3it /2. Contours from 0.0 to 9.6 x 10 -1+ with an interval of 
6.0 * 10 -5 . The labels are scaled by 10 6 . 
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Figure 53b. Contours of the phase, in degrees, of u of the first harmonic 

o, i. 

of the free stream oscillation of the flow past the parabola 
with a = 3tt /2 , contours from -280.0° to 60.0° with an interval of 
20 . 0 °. 
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1 mm 


Figure 54. Cross-section of the nose of the body used in the experiment of 
Kachanov, Kozlov, and Levchenko (1978). One inch on the figure 
is equal to 1 mm on the body; see the scale on the figure* 

This figure shows the body from the nose to a distance of 5mm 
back. The dashed line in the figure is the line y - o for both 
the upper and lower ellipses. The upper surface in this figure 
is that above which the measurements were made* 




Figure 55. Experimental results of Kachanov, Kozlov, and Levchenko (1978). 

Contours of u Q Q (in the notation of this paper) are shown. 

Contours are labeled with the speed in m /sec and the dashed 
curve labeled 6 is the position of the top of the boundary 
layer. Distances are in millimeters. 
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The steady and secondary streaming flow components, f^ Q and 
^2 o» as a f unct i° n of ti, the dimensionless distance 'from the 

boundary. The Strouhal number, a = tt/4. The solid curves are 
the results of numerical integrations and the dashed curve is the 
low frequency approximation for fl 



136 


1.0 -t 



- 0.6 - 
- 0.8 - 


-1.0 -J 

Figure 2. The 0(e), and and 0(e 3 ), ^21 and §2 1» 

amplitudes of the fundamental component of the oscillating flow 
as functions of n, the dimensionless distance from the bound- 
ary. The Strouhal number, a = tr/4. The solid curves are the 
results of numerical integrations and the dashed curves are the 
low frequency approximations for f^ ^ and g^ The low fre- 
quency approximation to is so close to the result of the 

numerical integration that the two curves are virtually indistin- 
guishable on this figure. 
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Figure 3. The 0(e 2 ) amplitudes of the first harmonic of the oscillating 
flow as functions of p, the dimensionless distance from the 
boundary. The Strouhal number, o = tt/4 . The solid curves are 
the results of numerical integrations and the dashed curve is the 
low frequency approximation to 2 ’ 


— v « KJ 


f J ^ x 50 ( low frequency approximation ) 


- 0.8 - 
- 1.0 — 

Figure 4 


The 0(e 3 ) amplitudes of the second harmonic of the oscillating 
flow as functions of n, the dimensionless distance from the 
boundary. The Strouhal number, a - tt/ 4. The solid curves are 
the results of numerical integrations and the dashed curve is the 
low frequency approximation to f’ 3. 
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Figure 5. The steady and secondary streaming flow components, and 

f l , as a function of ti, the dimensionless distance from the 
z , o 

boundary. The Strouhal number, o = 8ir. The solid curves are 
the results of numerical integrations and the dashed curve is the 
high frequency approximation to Q . 
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Figure 6. The 0(e) amplitudes of the fundamental component of the oscil- 
lating flow as a function of n, the dimensionless distance from 
the boundary. The Strouhal number, a = 8ir. The solid curves 
are the results of numerical integrations and the dashed curves 
are the high frequency approximations. The high frequency 

approximation to f' ■, is so close to the result of the numerical 

0,1 

integration that the two curves are virtually indistinguishable 
on this figure. 
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Figure 7. The O(e^) amplitudes of the first harmonic of the oscillating 
flow as a function of n, the dimensionless distance from the 
boundary. The Strouhal number, a = 8ir. The solid curves are 
the results of numerical integration and the dashed curves are 
the high frequency approximations. 
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Figure 8. The O(e^) amplitudes of the fundamental component of the 

oscillating flow as a function of ti, the dimensionless distance 
from the boundary. The Strouhal number, o = 8tt. The solid 
curves are the results of numerical integrations and the dashed 
curves are the high frequency approximations. 



v V 


Figure 9. The O(e^) amplitudes of the second harmonic of the oscillating 
flow as a function of n, the dimensionless distance from the 
boundary. The Strouhal number, a = 8it. The solid curves are 
the results of numerical integrations and the dashed curves are 
the high frequency approximations. 
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ABSTRACT 


A solution of the Navier- Stokes equations is given for an incompressi- 
ble stagnation point flow whose magnitude oscillates in time about a con- 
stant, nonzero, value (an unsteady Hieraenz flow). Analytic approximations 
to the solution in the low and high frequency limits are given and compared 
to the results of numerical integrations. The application of these results 
to one aspect of the boundary layer receptivity problem is also discussed. 
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1 . INTRODUCTION 


In this appendix we give a solution to the Navier-Stokes equations for 
an incompressible stagnation point flow whose magnitude oscillates in time 
about a constant, nonzero, value. Apart from the intrinsic interest of this 
problem, its solution is the first step in the solution of one aspect of the 
boundary layer receptivity problem, that is, the determination of the magni- 
tude and form of the disturbance introduced into the boundary layer on a 
body by a perturbation in the free stream. The solution of this problem 
would permit the calculation of the initial amplitudes of the Tollmien- 
Schlichting eigenmodes and continuum eigenfunctions in the boundary layer 
and give a rational foundation to transition prediction methods. 

We have recently given the solution to the boundary layer receptivity 
problem within the context of incompressible, linear stability theory for a 
parallel shear flow (Salwen, Kelly, and Grosch, 1980; Grosch & Salwen, 1980; 
Salwen and Grosch, 1981). There is, however, one aspect of the boundary 
layer receptivity problem to which our parallel flow solution is clearly not 
applicable. If we consider the flow near the forward stagnation point of a 
body, a linearized parallel flow theory cannot be valid because the flow is 
intrinsically nonlinear and nonparallel. Although we may be able to use the 
boundary layer equations away from the stagnation point, the full Navies 
Stokes equations must be used in the immediate vicinity of the stagnation 
point. Once a solution of the Navier-Stokes equations for the perturbed 
stagnation point flow has been found, it is possible to extend this solution 
away from the stagnation point and around the nose of the body. This is, in 
fact, what is normally done for the steady flow past an object. The Blasius 


A 
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series (Schlichting, 1979, pp. 168ff) is the extension around the nose of a 
blunt body of the Hiemenz solution (Schlichting, 1979, pp. 95ff) to the 
Navier-Stokes equations near a stagnation point. In section 6, we discuss 
the analogous extension for our oscillating flow solution. 

We were stimulated to do this work by the papers of Glauert (1956) and 
Stuart (1966). There are a number of other relevant studies in this area 
including those of Lighthill (1954) , Rott (1956) , Lin (1956) , Carrier and 
DiPrima (1957), Gibson (1957), Watson (1959), and Sarma (1964). Riley 
(1975) and Schlichting (1979, chapter 15) give comprehensive reviews of 
recent work in unsteady boundary layer theory. Lighthill (1978) has 
recently reviewed the current understanding of the phenomenon of acoustic 
streaming; i.e., the generation of a steady flow by the Reynold's stress 
due to an oscillating flow. 

2. EQUATIONS AND BOUNDARY CONDITIONS 

We consider the flow in the neighborhood of a stagnation point at 
(0,0) on a plane wall, with x the coordinate along the wall and y the 
coordinate perpendicular to the wall. We seek the solution (u(x,y,t), 
v(x,y,t)] of the two-dimensional Navier-Stokes equations which corresponds 
to the potential flow 

U = (U x/fc)(l + e cos w t ) (1) 

o 

V = - (U yA)(l + e cos tot) (2) 
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in the far field. Here, U and V are the x- and y-components of velocity 
of the potential flow, the velocity scale, i the length scale, £ 

the dimensionless amplitude of the oscillation, and to the frequency of 
oscillation. Defining a stream function 4>(x,y,t) by 


_ 3^ Sip 

cTy 5 V <Tx 


(2a, b) 


and substituting into the incompressible Navier- Stokes equations we find, as 
usual, that ^ is the solution of 


8 8iJ) 8 

'cTt + <Ty 5"x ~ 


V 2 4) = vV 2 (V 2 40 , 


(3) 


with v the kinematic viscosity and 

2 8 2 a 2 

V = + . 

8 x 2 8 y 2 


(4) 


We define dimensionless variables 


? = x/£, n = ( y/£ y R , T = U t/£, 

o o 


(5 a ,b , c) 


a Reynolds number. 


R = U l/v (6 a) 

o o 


150 



and a Strouhal number 

a = aslL/\j . (6b) 

o 

The use of separate velocity and length scales, and &, while 

correct, is somewhat arbitrary in that they appear only in the combination 

U A in the potential flow. The scales which are intrinsic to this problem 
o 

are the time scale of the base flow, T = i /V , the frequency, <o, and 

o o 

the kinematic viscosity, V. We can define velocity and length scales in 
terms of T q and V, thereby setting the Reynolds number equal to unity. 

In any case, the dimensionless parameters which appear in the equations and 
boundary conditions are the oscillation amplitude, £, and the Strouhal 
number, a, and the results will be the same. We have chosen to use 
independent velocity and length scales, U q and i, because we intend to 
use the solution presented here as the basis for constructing a solution to 
the problem of an oscillating flow past a blunt object with length scale A. 

If we set, in analogy with the Hiemenz solution ( Schlicht ing, 1979, pp. 
95 f f ) , 


= (U £// R ) K F(n ,t) , 
o o 


(7) 


then equation (3) reduces to 


*— [3ii_ * (iLl) 2 - F *1 - ill] - o 

3n 3x3n 3h 3n 2 3h 3 


( 8 ) 
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which may be integrated to yield 


3 2 F 
3x3ti 


3n 


3 2 F 3 3 F 

3n 2 3 ti 3 


G(5 ,T ) . 


(9) 


Because the left-hand side of equation (9) is independent of ?, it is 
clear that the "constant of integration", G, only depends on the 
dimensionless time, x. In order that the x and t dependence of the 
pressure in the far field agree with that of the potential flow, we must 
have 


(U /A) 2 G = + = (U /S,)(-eia sm wx) + (U /£) 2 

o ~ ~ o ' o 

x 3q x 3 x 


x (1 + £ cos utf . 


( 10 ) 


Therefore, F(h ,T ) is the solution of 


3 2 F .3F 2 3 2 F 9 3 F 

+ ( — ) - F 

3x3n 3h 3n 2 9n 3 


= 1 + e(2 cos ax - a sin ax) 


+ j e 2 (1 + cos 2ax ) 


(ID 


with the boundary conditions 


F (0 ,x ) = 0, 


( 


3F 

5rf'n =0 


0, 


(12a) 


(12b) 


j 
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•=— + 1 + cos ot as n 

an 


(12c) 


3. METHOD OF SOLUTION 


To solve equation (11), we will expand F(ti,t) in a Fourier series in 
the dimensionless time, t. The coefficients of the expansion are functions 
of n and the parameter e and each of these will be expanded in a power 
series in e. It is easy to see, from the form of (11), that the T- 
independent term in the Fourier series contains only even powers of e. 
Therefore, we look for a solution of the form 


F(u,t) = £ (e m f (n) + J e [f (n) cos kar 

2m, o , , 2m, k 

m=o k=l 

+ g 2m k (h) sin kcrr]}. (13) 


To find the equations obeyed by the f 2m k and ^ we substitute 

(13) into (11), collect the coefficients of like terms in the Fourier 
series, and set the coefficients of the successive powers of e equal to 
zero. We find, first, that f is the solution of 


d 3 f d 2 f df 2 

+ f - ( — °iSj) +1 = 0 (14) 

dr|3 °»° ^2 dri 


with 
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f (0) - f' (0) = 0, f' + 1 as n + «. 

0,0 0,0 0,0 


(15a, b,c) 


As was expected , f , the t- independent , e- independent term in the 

0,0 

series, is the Hiemenz solution for the steady stagnation point flow. 
Next we define the operator L by 


L = — + f — - 2f ' — + f" , 

dn 3 °’° dn 2 °’° dn 


(16) 


with primes denoting differentiation with respect to n. Then it can be 
shown that the equations and boundary conditions for the ^ an< * 

{ g 2m,k* 

are, for k > 1, 


Lf - kff g 1 = P , 

2m, k 6 2m, k 2m, k 


Lg + kO f 1 = Q , 

s 2m,k 2m, k v 2m,k’ 


(17b) 


f_ . (0) = g_ . (0) = f' (0) = g' , (0) = 0, 
2m, k 2m, k 2m, k 2m, k 


(17c) 


f' ->-6 6 , g' 0 as h 00 , 

2m, k m,o k,l 2m, k 


( 1 7 d ) 


and, for m > o and k = 0, 


L f. = R„ 

2m, o 2m, o 


(18a) 
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f, (0) = 

2m, o 


f', <0) = 

2m, o 


0 , 


(18b) 


f' -*-0 as ti (18c) 

2m, o 

The { P { Q are linear combinations of products of the 

2m, k 2ra,k 

{ f_ } , { g„ } and their derivatives, where s < k and 0 <- r m or 

2r,s 2r,s 

s = k and r < m, and are given in Section 7. Therefore, these equations 
can be solved sequentially. The equations can be integrated numerically 
quite easily, although care must be taken to control roundoff errors. The 
results of these integrations are given and discussed in Section 5. 

4. LOW AND HIGH FREQUENCY APPROXIMATIONS 
4.1 The Low Frequency Limit 

As o 0, it is expected that the solution will approach a quasi- 
steady solution. It is straightforward to show that the quasi- steady 
solution. 


F(n ,t ) = (l + e cos ai) 1/2 f ([ 1 + e cos at l 1 /2 n) , 

0,0 


(19) 


satisfies equation (11) and the boundary conditions (12) to 0(o). 

It is also easy to show that this quasi- steady solution is consistent 
with the expansion given in equation (13). If F, as given in (19), is 
expanded in a Taylor series in e , we find that 


F = l 


(e cos ot ) 


n=o 2 n! 
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00 oo 2nH-k S . (h) 

1 <2 - 5 ko > £ <f> — ^ - cos k0T > 

k=o m=o m ! (m + k) ! 


( 20 ) 


with 


s (n ) = f (n ) , n = o 

n o,o 


= n 


n f (n) (n ) + l (-1) 


n- 1 n- r- 1 


o ,o 


r=o 


(2n - 5 - 2)1 
2 n ' r ' 1 r!(a-r)> 


[o - r(r 


n > 1, 


( 21 ) 


and 


, v d f 

f (r) = o»o (22) 

°’° " dn r 


Equating the coefficients of e in (21) to the corresponding 
coefficients of £ in (13) we find that, to order o. 


C 2m,0 2 4m (tn!) 2 S2m ’ 


(2 3a) 


L 2 m ,k 


(23b) 


’2m,k 


= 0 


(23c) 
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Taking the a -*• 0 limit in the differential equations (17) and (18), for 


the {f } and {g }, and substituting the low frequency approxima- 

4in j k zm, k 

tions (23), it can be shown by induction, after some lengthy but straight- 
forward calculations, that these approximations satisfy the dif ferentional 
equations to 0(a). Therefore, the quasi- steady solution is the same as the 
expansion in the low frequency limit. 


4.2 The High Frequency Limit 


In the high frequency limit, a 00 , it would be natural to look for 

approximate solutions for the {f ,} and { g„ ,} in the form of a power 

2m, k 2m, k 

series expansion in a - *- . However, it is clear from the form of the 
differential equations, that this expansion would be non-uniform because the 
highest derivatives would be multiplied by the small parameter o - ^ and 
thus would vanish as a °° . 

We therefore rescale the equations, defining a new independent 
variable, an inner variable, 


z = (a/2) 1 ^ 2 rt = (w/2v) 1/2 y. 


(24) 


■j / O 

with the length scale (2v/w) , that of Stokes's second problem 

(Schlichting, 1979; pg . 75 and Chapter 11). In this limit we assume that 
there is an inner expansion of the form, 


= f (n) 

o ,o 


I 

rrF 1 


2m "2 
e f. 


2m, o 


(z) 


oo oo 

+ I I e 2m k [ cos kax + g (z) sin kax] (25) 

zm, K zm, iC 

m=o k=l 
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It is easy to see that the differential equations for the {f„ , } 

zm, k 

and {§2m are j ust ec l uat i ons (17) and (18) with a transformation of 
variables from t| to z. It can be inferred from the form of the equations 
that fj^ k and g2 m ^ are 0((2/a)^ m+ ^ ^^) in the high frequency limit. 
We then solve these equations, retaining only the highest order terms in an 
expansion in powers of (2/a). 

We find 

1/2 

f Q ^(z) = (2/a) [z + y e 2 (cos z - sin z) - — ], (26a) 

g Q ^(z) = (2/a)l /2 [~ e 2 (cos z + sin z) - -i], (26b) 

a Stokes shear layer flow caused by the 0(e) part of the far flow field. 

In the high frequency limit it is decoupled from the 0( e°) steady outer 

flow. This is, of course, a familiar result in time dependent boundary 

layer theory and has been derived and discussed by Carrier and DiPrima 

(1957), Stuart (1966), and Riley (1975), among others. 

If we let 6 and 6 , be the boundary layer thickness of the 

0,0 o,l 

steady flow and the Stokes shear layer flow, respectively, then, for 
a » 1, 

6 /6 « 2 (2/a) 1/2 « 1. (27) 

0,1 0,0 

Next, we solve the equations for the 0(e 2 ) oscillatory flow and find 

that 
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f „(z) = (2/a ) 3/, 2{!^ [ e ^ Z (sin z + cos /T z) - l] 
° jZ 16 


1 -z . , 

+ 2 j- ze sin zj 


(2 8a) 


g _ (z) = (2 /a) 3/2 {_?. [e ^ (sin /T z + cos /T z) - l] 
°* 16 


1 _z l 

ze cos z) 

4 


(28b) 


Again we have a Stokes shear layer, decoupled from the steady flow and, for 
large a , confined to a thin layer imbedded within the steady boundary 
layer . 

Proceeding next to the 0(e 2 ) steady streaming flow component, we find 
that the most general solution which satisfies the boundary conditions at 


z = o, is 


f 0 = (2/a ) 3/2 (i_) [ 13 - 6z - e - 4e (3 cos z + 2 sin z) 

2 ’° 16 


- 4ze Z sin z] + Dz 2 , 


(29) 


where D is an arbitrary constant. This secondary steady streaming flow is 
identical to that found by Stuart (1966) using the boundary layer equations. 
Stuart's small parameter, a, is the reciprocal of a and Stuart's expan- 
sion is in powers of a, while ours is in powers of (2/a). Therefore, 
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A , # . f 

f , with D - 0, is equal to 1/2 of the x given by Stuart (see Stuart, 
2,o s 

1966; eq. 2.11). Our solution and Stuart's satisfy the boundary conditions 

at the wall but do not satisfy the outer boundary condition because 


df 


2,o 


dn 


- 


(2/a) + a Dn as n -*■ 


(30) 


and does not vanish, as required by the far field boundary condition, for 

any value of D. This is explained by the fact (pointed out by Stuart) 

that -there is an outer boundary layer, thick compared to the Stokes layer 

thickness, 6 , but thin compared to the scale of the body, within which 

0,1 

this secondary steady streaming flow decays to zero. 

Since the a °° limit is non-uniform and z is an inner variable, 

A A 

all of the f 1 s and g' s are inner solutions and could be expected to 

A. 

require matching to appropriate outer solutions. The functions f , 

o, 1 

A A A , 

g , f , and g , however, are also outer solutions and, unlike 
o,l o,2 o,2 

A. 

f , do not require matching. 

2,o 


Outside of the Stokes layer, the interaction between the secondary 

A A 

streaming flow, f' , and the oscillatory potential flow, f' , is 
2,o o, 1 

unimportant. In the case studied by Stuart, the dominant non-linear 

A. 

interaction is that between f^ and itself. Since the velocity of this 

flow at the top of the Stokes layer is 0(U o /a), the corresponding Reynolds 

number is R = (U /a)(£/v) = R /a. From this, Stuart concluded that the 
so o 

thickness of the outer layer is 0(^R = 0(£R - 1/2 a 1 ^ 2 ), which is 

so 

larger than that of Stokes layer by a factor of the order of a. 
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For the case under consideration here,, the situation is completely 


changed by the existence of the large steady flow, f^ q with boundary 

layer thickness 6 = 0 (a 1/<2 6 ). The dominant non-linear interaction 

o,o o,I 

of f 0 in the outer layer is with f' and, consequently, it falls off in a 
2 , o o,o n J ’ 

distance of the order of 6 (instead of o l/2 <5 ). The appropriate 

o,o o,o 

outer variable is therefore ri = y/6 

0,0 


If we express f , 
0,1 

and 

g 1 in terms 
0,1 

of n. 

substitute them into 

the differential 

equation 

for 

f„ , and let 
2,o 

a *■ °°, 

it can be seen that the 

outer solution. 

f ? is 
2,o 

the 

solution of 





) 

o 


3 2 

(- + f 2f ' — + f" ) f = 0, 

dn 3 °»° d n 2 °’° dn °»° 2 >° 


(31) 


with the outer boundary condition 


f'» o as h oo. 

2,o 


(32) 


Using the asymptotic expansion for f Q q . 


f ~ n + A, 
0,0 


(33) 


with A a constant, it can be shown that the general solution of (31) is 


f 0 n = C f* (n) + C.h, (n) + C„ h_( 
2,o o o,o 1 1 y 22 


n) 


(34) 


where, as n 


~ (i + e 2 ) / 


e~ Y /2 dr - He '* 2/2 


(35) 
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h' 


2 




8 


+ 



with 


(36) 


8 = n + A. 


(37) 


The outer boundary condition requires that C 2 be set equal to zero. 

This leaves three arbitrary constants, C , C , and D, in the inner and 

o 1 

outer solutions so we can match the inner and outer solutions and their 
first two derivatives. Matching the inner limit of the outer solution and 
the outer limit of the inner solution shows that 


C q = O(0 _1 ) , Cj = 0(cr- 3/2 ), and D = 0(o" 2 ) . 


(38a, b, c) 


The inner solution was obtained by expanding in powers of a - -*-^ 2 and 
retaining only the lowest order terms. Consistent use of this approximation 
requires that we set and D to zero and match the first derivative. 

We find that 


C = -3/(4B0) . 
o 


(39) 


where B — f" (0). The composite solution for f is 

0,0 2,o 


^2 0 ~ (2/<j) 3 ^ 2 (-j^-) [ 13 - e - 4e (3 cos z + sin z) 
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(40) 


- 4ze z sin z] — f' (q), 

4Ba 0,0 

which satisfies the boundary conditions at z equal to zero as well as 

z -*■ °°. The thickness of the inner layer of the steady streaming flow is 

6 i , while that of the outer layer is <5 which is much larger than 

0,1 0,0 

6 . for a 00 . 

o,l 

We have also calculated the high frequency approximations for the 
0(e 3 ) functions, for f Q g Q f 2 and These are given in 

Section 8. The components of the amplitude of the 0(e 3 ) portion of the 

second harmonic f „ and g - are driven by the interaction of the 0(e) 

fundamental (f , ,g . )and the 0(e 2 ) first harmonic (f OJ g _ ) . The 
0,1 o,i o , 2 o , 2 

inner expansions for f „ and g „, given in Section 8, are also outer 

O j J O j J 

expansions . 

This is not true for f^ ^ and g^ the components of the amplitude 

of the 0(e 3 ) part of the fundamental. They are excited by the interaction 

of (f ., g .. ) with (f_ ) as well as (f , ,g ,) with (f „,g „). 
o,l o,l 2,o o,l o,l o,2 o,2 

The inner expansion of the in-phase component, f , satisfies the outer 

boundary condition, but inner expansion for the out-of-phase component, 

g , does not. The outer expansion for g„ 1 is the same as that for 

2,1 

f 2 Q . The composite expansion for g 2 ^ is found in the same way as that 
for f 2 Q and is given in Section 8, 

5 . RESULTS 

We have numerically integrated the differential equations for the f's 
and g's with 2m + k < 3 over a wide range of values of a. Some typical 
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results at a moderately low frequency, O = ir /4 , are shown in Figures 
1-4, where the solid curves show the results of the numerical integrations 
and the dashed curves are obtained from low frequency approximations for the 
f's given in section 4.1. The g's are identically zero in the low fre- 
quency approximation. 

It can be seen from the results shown in these figures that at low 

frequencies there is no Stokes layer; the boundary layer thickness of the 

various components are generally equal to or greater than that of the steady 

flow component, f' . The steady streaming component, f' , is quite 

0,0 2 ,o 

small compared to the mean flow, even for e = 1. 

At a=o the g's, the out-of-phase components, are identically 

zero. At small, but non-zero, frequencies the low order (in e) in-phase 

components, the f's are larger than the g's. As the order increases, 

however, the magnitudes of the f's and g's tend to equalize and 

decrease. 

It is also apparent from an examination of Figures 1-4 that the low 

frequency approximations are reasonably accurate even at o = tt/ 4. In fact, 

the difference between the numerical solution and the approximation for 

f is so small that it is not apparent in Figure 2. Taking into account 

o,l 

the fact that the higher order terms, which have the largest deviations from 
the low frequency approximations, have very small magnitudes, it is clear 
that the low frequency approximations, equation (19), is reasonably accu- 
rate, even for Strouhal numbers as large as Tf/4. 

Figures 5-9 show the f's and g's for 2m + k < 3 at a moderately 
high Strouhal number, o = 8n , as obtained from the numerical integrations 
(solid curves) and the asymptotic approximations (dashed curves). The top 
of the Stokes layer is at z = 4.6 and this corresponds to n “ 6.5 a . 
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Therefore, the top of the Stokes layer is at n “ 1.3 for a = 8 tt. The 
secondary streaming flow (see Fig. 5) extends far beyond the Stokes layer. 
The variation of the fundamental component (see Fig. 6) is confined to the 
Stokes layer. The amplitude functions of the higher harmonics tend, like 

f 2 q , to extend beyond the Stokes layer because they are also driven by 

Reynolds stresses due to f 

J 0,0 

The secondary steady streaming flow, f' , is considerably smaller 

/ , o 

at high frequencies than at low frequencies. In addition, the nature of the 
secondary flow changes as a increases; at a small a, the net secondary 
flow is positive while, for large o, it is negative. As discussed by 
Stuart (1963), it is known from experiments that this effect occurs in the 
steady streaming flow generated by a circular cylinder oscillating along a 
diameter. 

3F 

The tangential velocity component, u, is proportional to -t^, that 

is to the ( f 2 m §2m k^ ‘ -*- n ^8 ure ^ we present some of the results of 

the numerical integrations; a plot of the variation of max If' I and 

zm , k 

max |g' , as functions of a. The maximum of f' is 1.0 for all a 

1 2m , k 1 o,o 

and the absolute maximum of f' , is about 1.069 at a « 17 ; the maxima of 

o,l 

f' and f' , are not plotted in Figure 10. We can conclude, from the 
0,0 o,l 

results shown in Figure 10, that the high frequency estimate, { | f^ ^ | , 

ii ii n f -2m-k+l) . . , 

l g 2m,kl^ = °^ CT 5 1S 9 ul te good. 

6. DISCUSSION 

The applied far-field flow consist of a time-independent mean flow and 
a fundamental with a coswt time-variation. Non-linear interactions result 
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in the generation of components at all multiples of the fundamental fre- 
quency and modifications of the mean flow and fundamental components, 
though, at moderate £ ( 1), the fundamental and first harmonic tend to 
dominate. 

At low frequencies, our results are well represented by a quasi- steady 
approximation (19), which has the same form as the steady stagnation- point 
flow except that the amplitude and scale vary with time. Somewhat surpris- 
ingly, this approximation is quite good for a dimensionless frequency, o, 
as large as tt /4 . 

At high frequency, the viscous boundary layer corresponding to the 

oscillating components is largely, but not completely, confined to a Stokes 

layer of thickness (2 /o) 1 ^ times that of the steady boundary layer. For 

these large values of o , the inner asymptotic approximations are solutions 

of differential equations which are independent of f and are, there- 

o, o 

fore, decoupled from the mean flow. For a number of components, the inner 

expansions are also the correct outer expansions; these components are, 

therefore, totally decoupled from the mean flow to lowest order in a and 

are, in fact, identical to the corresponding solutions found by Stuart 

(1966) for a purely oscillatory flow past a body. 

Not all of the components are decoupled from the mean flow in the high 

frequency limit. The high frequency expansion is a non-uniform asymptotic 

expansion and thus, for certain components, in particular the steady second 

order flow, f' , an outer expansion, matched to the inner expansion, is 
z , o 

needed for a uniformly valid approximation. The results given in Section 
4.2 show that such composite expansions can be found for a >> 1; they 



satisfy, to lowest order in a, the differential equations and all of the 
boundary conditions. We believe, although we have not carried out the cal- 
culations, that this is also true for an expansion to any order in a. The 
results of the numerical integrations, shown in Figures 5-9 are consistent 
with this belief. We conclude that, both for this flow and for the oscilla- 
tory flow studied by Stuart (1966), there are second, and higher, order 
steady streaming flow components that do not extend to infinity but are 
confined to a layer adjacent to the solid boundary. This layer is much 
thicker than the Stokes layer but much thinner than the length scale of the 
body . 

There are other aspects of Stuart’s solution which can be compared with 

ours. He used three parameters, a, B, and R in his solution, which 

s 

are related to our Strouhal number, o, and Reynolds number, R , by: 

o 

a = u /(goJI ) = a- 1 , 
o 

B + v /(coil 2 ) = 1 / (oR q ) , 

R = u 2 /(rnv ) = r /a . 
so o 

Stuart used the boundary layer equations and an expansion for ct -*■ o; we 

have used the Navier-St okes equations and, in addition to the numerical 

solutions, an expansion for the high frequency limit, o 00 . It is clear 

that Stuart's expansion, and ours, are high frequency approximations. 

Stuart assumed that B was small and R large in order to justify the use 

s 

of the boundary layer equations. In the high frequency limit a 00 , and, 

with R fixed, a, B, and R are all small. We do not need to make 
o s 

any assumptions concerning the magnitude of R q or R g because we used the 
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Navier-Stokes equations. We can determine the f's and g's by numeri- 
cally integrating the ordinary differential equations; however, the quasi- 
steady solution for low frequencies and the high frequency expansion are 
useful analytic approximations. 

We have not attempted to determine the radius of convergence, in e, 

of the series for F(q,T,e), equation (13). However, some observations can 

be made. We have found that f' , and g„ , are 0(a ^ f or 

2m , k ° 2m , k 

large a. Therefore, max { | f' |,|g' |}, the amplitude functions for 

Za\\ j K zm j K 

-2m-k+lv _ . . _ , . s 2m+k 

u, are 0 (.0 ;. Thus the series for u is in terms of (£/a; , 

for large a, and this suggests that it converges for e/a < 1 and 

converges rapidly for e < 1. On the basis of the results presented in 

figure (10) we conclude that the high frequency bound on max {If' I, 

2m, k 

|g' 2 m is valid at all frequencies and, therefore, that the series 

converges rapidly for e < 1. Of course the convergence will be slower at 
low frequencies but we have shown that, for a ~ tt/ 4, the quasi-steady 
solution is an accurate approximation. 

In the introduction we suggested that the solution to the problem of 
oscillating stagnation point flow was the first step in the solution of one 
aspect of the boundary layer receptivity problem. A few years ago, Morkovin 
(1978) reviewed the rather rudimentary state of knowledge of the dynamics of 
boundary layer receptivity. He identified four general classes of free 
stream disturbances which might generate Tollmien-Schlichting waves in the 
boundary layer. These are: vorticity fluctuations, sound, entropy 

disturbances, and unsteadiness in the mean flow. In the incompressible 
limit, there can be only vorticity fluctuations and unsteadiness. Morkovin 
argues, however, that a temporally oscillating incompressible flow is a 



reasonably accurate approximation for a sound wave impinging normally on the 
nose of a body if the wavelength of the sound wave is much greater than the 
radius of curvature of the nose. 

The solution of the Navier-Stokes equations given here is the solution, 

in the neighborhood of the stagnation point, of the receptivity problem for 

a simple unsteady flow. In order to interpret this solution in terms of a 

stability model it is necessary to extend it around the nose of the body 

into a region where the flow is, at least, quas i- paral le 1. If we can assume 

that the tangential component of the potential flow at the edge of the 

boundary is of the form U = U H(5) (1 + e cosat), with H(£) having a power 

o 

series expansion in £ , the distance along the body from the stagnation 
point, we can use an unsteady variant of the Blasius series to solve this 
problem. For the first term we use the full two-dimensional Navier-Stokes 
equations and the solution is that given here. For the subsequent terms, we 
use the boundary layer equations. We have carried out this calculation 
procedure up to and including the terms in C 1 * and we are now applying it 
to bodies which are elliptic or parabolic cylinders. We hope to report the 
results of these calculations at a later date. 

This work was supported, in part, by grant NAG- 1-96, from the National 
Aeronautics and Space Administration. 
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7. DERIVATION OF P. . AND Q„ . MATRICES 

2m, k 2m, k 

The expressions for P„ , , Q„ , , and R. in (17a), (17b), and (18a) 

r 2m, k 2m, k 2m, o 

are simplified by use of the notation 


M(f,g) = f"g - 2f ' g ' + fg". 


(7-1) 


In terras of M, 


P ok - -i t85 kl + 2S M + (1 - V X [ M(f o,r> f o,k-r> 

r=l 


‘ M(g o,r' g o,k-r)l } f ° r kil> 


(7-2) 


1 | ra-1 m k-1 

P 2m,k “ 4 [ 4 ^ M(f 2(m-s) ,o ,f 2s,o ) + (l " 6 kl } ^ ^ 

5 i s ’ s =o r=l 


s=o 


^^2s ,r’ ^2(m-s) ,k-r^ ^ g 2s ,r ,g 2(m-s) , k-r^ 


m-1 m-l-s 

+ 2 l l [m( f 


, 2s , k+m-s-s ' * ^2s ' ,m-s-s ' ^ ^ g 2s , k+m-s-s 1 5 

s=o s =o 


g 2s,m-s-s'^ f ° r m - 1 * k - X » 


(7-3) 


k-1 


Q ok = 1 2a6 kl - (1 - 5 - } ^ M(f ~ - 


kl‘ 


) M(f ,g , ) 

, o,r o,k-r 
r=l 


for k — 1, (7-4) 
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kl 



and 


Here, 


Q 2m,k 2 2 l M ( f 2(m-s) ,g 2s,k 

’ s=o 


m 


S=i 


m-1 m-l-s 


W f 2s,r’ g 2(m-s) ,k-r^ + ^ ^ f 2s, 

5 s=o s =o 


X g 2s 1 ,k+m-s-s' ’ f 2s ' ,k+ra-s-s 1 ’ S 2s 1 ,m-s 


for m-1, k - 1, 


2m 


m-1 

= - 7 2i , + 2(1 - « .) I M(f ,f f \ 
, o 4 ml ml , 2s, o 2(m-s) 

5 s— I 


m-1 m-l-s 

l ,1 t M ( f 2s,m-s-s , ’ f 2s , ,m-s-s'^ + M ^ { 


s=o s =o 


X 


g 2s ' ,m-s _ s ' 


)] 


for 


> 

m - 


1 . 



is the Kronecker 6 symbol. 


k-1 

I 

r=l 


m-s-s ' 



(7-5) 


) 

,o 


2s , m-s-s ' , 


(7-6) 
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8. HIGH FREQUENCY APPROXIMATIONS TO THE 0(e 3 ) FUNCTIONS 


• A 

The approximations to f _ and g _ for a + 00 are: 

o , j 0,0 


f Q 3 ( z ) = (2/a) {[9 - 8/2 + 3 (2 + l)] + e" Z [(/T 


2 2 
-6 + 4 z - 4z ) cos z + (6 - /2" - 4z ) sin z] 


+ e Z ( 4 + 2/r - 8z) cos /2z - (4 + 2/2) s in /2z] 


- / 3(/"2 + l) e Z ( cos /3z - sin /*3zJ + ( 5/T - 7 ) 


x e^ [ cos (/IT + l) z - sin ( 2 + l)z]} /128, 


(z) = (2/a) {[ 9 - 8/2 + /3 (/T + l)] + e 2 [(/T 


2 2 
-6 4z J cos z + (/~2 - 6 + 4z - 4z ) sin z] 


+ e Z [ ( 4 + 2/2) cos /2z + (4 + 2/7 - 8z) s in /2z] 


- /3[/T + l) e ( cos /3z + sin /3z) + ( 5 /T - 7) 


x e ^ ^ Z [ cos (/T + l) z + sin (/T + l)z]}/128. 
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The approximation to f„ 1 for 0 -*■ 00 is: 


5/2 


f (z) = (2/o) {( 349-335/7) /1920 + e“ Z [ ( 60/2-213) /9 60 

Z , 1 


+ ( (5z/64) (z-l)/64) cos z + ((93 + 15/2) /240- lOz 


x (z + 1)/64J sin z ] + e” [ (7/2" - 2 + 4z) cos /2 


- (7/2" + 2) sin/2z]/64 + (13/160)e“ 2z + e _( ^ 2 * +1)z 


[(/T-5) cos(/T- 1)z - (/7 + 5) sin (/T - l)z]/384 


3 z 

+ e cos z/320} 


(8-3) 


The composite approximation to g for 0 “ is: 

2,1 


5/2 


g (z) = (2/0) { - ( 1513 + 125/2) /1920 + e~ Z [(558- 45/2) /9 60 

Z , 1 


+ ( (3z (6-z)/64) cos z + ((177 + 45/1) /960 + 


x z(31 + 9z)/64) sin z] + e (7/2" + 2 ) cos /2 


-2z r 


+ ( 7 / 2 "- 2 + 4z) sin /2z] /64 + e [(13/80 + z/16)] 
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z[ (5/2) cos(/T - l)z - (5 - /2) 


with 


-(/2 + 1 ) 

+ e 


x sin (/2 - 1) z] /384 + e Z s in z/320} 

2 

+ ( 2 /a) 3f ' (n)/(8B), (8-4) 

0,0 

B = f" (0). (8-5) 

0,0 
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Figure 1. The steady and secondary streaming flow components, f^ and 

f 2 j0 > as a function of the dimensionless distance *?rom the 

boundary. The Strouhal number, o = tt/ 4. The solid curves are 
the results of numerical integrations and the dashed curve is the 
low frequency approximation for f£ Q . 
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1.0 



Figure 2. The 0(e), f^, and g^, and 0(e 3 ), f£ t and j, 

amplitudes of the fundamental component of the oscillating flow 
as functions of h, the dimensionless distance from the bound- 
ary. The Strouhal number, or * ir/4. The solid curves are the 
results of numerical integrations and the dashed curves are the 
low frequency approximations for f^ ^ and g^ ^ . The low fre- 
quency approximation to is so close to fhe result of the 

numerical integration that the two curves are virtually indistin- 
guishable on this figure. 
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Figure 3. The 0(e 2 ) amplitudes of the first harmonic of the oscillating 
flow as functions of n, the dimensionless distance from the 
boundary. The Strouhal number, a = ir/4. The solid curves are 
the results of numerical integrations and the dashed curve is the 
low frequency approximation to 2’ 


i 
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Figure 4. The 0(e 3 ) amplitudes of the second harmonic of the oscillating 
flow as functions of n, the dimensionless distance from the 
boundary. The Strouhal number, a = tt/ 4. The solid curves are 
the results of numerical integrations and the dashed curve is the 
low frequency approximation to f' o. 
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Figure 5. The steady and secondary streaming flow components, and 

f I , as a function of n. the dimensionless distance from the 
z,o’ 

boundary. The Strouhal number, a = 8ir, The solid curves are 
the results of numerical integrations and the dashed curve is the 
high frequency approximation to f£ Q . 
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Figure 6. The 0(e) amplitudes of the fundamental component of the oscil- 
lating flow as a function of n, the dimensionless distance from 
the boundary. The Strouhal number, a = 8ir. The solid curves 
are the results of numerical integrations and the dashed curves 
are the high frequency approximations. The high frequency 
approximation to f^ ^ is so close to the result of the numerical 

integration that the two curves are virtually indistinguishable 
on this figure. 
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Figure 7. The 0(e 2 ) amplitudes of the first harmonic of the oscillating 
flow as a function of n, the dimensionless distance from the 
boundary. The Strouhal number, o = 8ir. The solid curves are 
the results of numerical integration and the dashed curves are 
the high frequency approximations. 
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Figure 8. The 0(e 3 ) amplitudes of the fundamental component of the 

oscillating flow as a function of n, the dimensionless distance 
from the boundary. The Strouhal number, a = 8 if. The solid 
curves are the results of numerical integrations and the dashed 
curves are the high frequency approximations. 
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Figure 9. The 0(e 3 ) amplitudes of the second harmonic of the oscillating 
flow as a function of n, the dimensionless distance from the 
boundary. The Strouhal number, a - 8ir. The solid curves are 
the results of numerical integrations and the dashed curves are 
the high frequency approximations. 
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Figure 10b. Variation of max | CT 1-2m_k f^ m fc | and max| <j l ~ 2m " k g^ m k | with 
the Strouhal number, a. 
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